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ON THE MOTION OF PARTICLES IN GENERAL 
RELATIVITY THEORY 


A. EINSTEIN and L. INFELD 


1. Introduction. The gravitational field manifests itself in the motion of 
bodies. Therefore the problem of determining the motion of such bodies from 
the field equations alone is of fundamental importance. This problem was 
solved for the first time some ten years ago and the equations of motion for two 
particles were then deduced [1]. A more general and simplified version of this 
problem was given shortly thereafter [2]. 

Mr. Lewison pointed out to us, that from our approximation procedure, it 
does not follow that the field equations can be solved up to an arbitrarily high 
approximation. This is indeed true. We believe that the present work not 
only removes this difficulty, but that it gives a new and deeper insight into the 
problem of motion. From the logical. point of view the present theory is 
considerably simpler and clearer than the old one. But as always, we must 
pay for these logical simplifications by prolonging the chain of technical 
argument. 

The subject matter is presented here from the beginning and the knowledge 
of previous work is not assumed. To facilitate the reading for those who have 
studied the previous papers we use here essentially the same notation as before. 

Let us start with some general remarks. 

All attempts to represent matter by an energy-momentum tensor are un- 
satisfactory and we wish to free our theory from any particular choice of such 
a tensor. Therefore we shall deal here only with gravitational equations in 
empty space, and matter will be represented by singularities of the gravitational 
field. 

In Newtonian mechanics, particles are represented as singularities of a scalar 
field y, which satisfies Laplace’s equation everywhere outside the singu- 
larities. Because the classical equation is linear, the field can be decomposed 
into partial fields, each part due to a single particle. Each particle is in a field 
due to all other particles. The theory is completed by the equation of motion, 
that is by putting the acceleration equal to the negative gradient of the field, 
the proportionality factor being a universal constant. Thus classical physics 
postulates the equations of motion independently of the field laws. The masses 
of the sources of the field are assumed to be independent of time. The laws of 
motion are supposed to be valid in an inertial system. Therefore space-time 
appears as an independent physical entity. The conceptual weakness of such 
a space-time background in the classical theory was already recognized by 
Newton. 
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If we compare this state of affairs with that in general relativity theory, in 
its original formulation, we see striking similarities and differences. Laplace’s 
equation 

Ag = 0 


is replaced by the gravitational equation 
Ri = 0, 


which, however, unlike the classical equation, satisfies the general relativity 
principle. The classical principle of inertia becomes in relativity theory the 
principle of the geodesic line valid for a particle with infinitely small mass. 
True enough, the difficulty with the inertial system disappears in relativity 
theory, as does the independent physical reality of space-time. Yet the equa- 
tions of motion still appear independently of the field equations. 

Our aim is to investigate to what extent the field equations alone contain the 
equations of motion of particles; also to develop a method that will allow us 
to find these equations of motion up to an arbitrary approximation. 

Let us start with a simple remark: a linear law always means that the motion 
of singularities is arbitrary. If to a world-line of a singularity with mass m, 
there belongs a field Fi.) and if to a world-line of a singularity with mass m2 
there belongs a field Fy), then the superposition of these two fields, that is 
F)+ Fq) is alsoa solution of the linear field equations. In such a solution the 
same two world-lines would appear together that before appeared singly. 
Therefore the field with its linear laws cannot imply any interaction between 
the singularities. Thus only non-linear field equations can provide us with 
equations of motion because only non-linearity can express the interaction 
between singularities. 

But the argument cannot be reversed. Non-linearity is necessary but not 
sufficient for the equations of motion to follow from the field equations. 

The reason why the gravitational field equations do provide us with equa- 
tions of motion lies not in their non-linear character alone, but also in the fact 
that these equations are not independent from each other Indeed, among the 
ten components four are free, this being due to the freedom of choice in the co- 
ordinate system. The ten equations are valid, so to speak, only for six effective 
functions. They would be inconsistent were it not for the four (Bianchi) 
identities that they satisfy. This must be so for every relativistic system of 
equations derived from a variational principle. These identities are (besides 
the non-linearity) responsible for the equations of motion being determined 
by the field equations. 

The ideas leading to the equations of motion are not easy and are mutually 
interwoven. 

One of the essential ideas in this paper is the treatment of gravitational 
equations by a “new approximation method.” In it we treat space and time 
differently. We regard the changes of the field in time as small compared with 
those in space. Only then do we arrive at a consistent, manageable set of 
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equations that can be solved step by step. This idea is not new and was con- 
tained in the previous papers. 

The other important idea is the deduction of the equations of motion, which 
are ordinary differential equations, from the field equations which are partial 
differential equations. This idea, treated here differently than in the previous 
papers, leads to the use of surface integrals taken around the singularities of 
the field. These surface integrals will depend only on the motion of the 
singularities and not on the shape of the surface. 

These and other ideas will be treated in detail in this paper. To make them 
clear we have decided to delegate all the more tedious calculations to the 
Appendices. (If we refer, for example, to A.4, this means the Appendix be- 
longing to Sec. 4.) But even so, many straightforward but long calculations 
had to be omitted. This is especially true for the calculations that lead beyond 
Newtonian motion. We included here a short section on this subject, just for 
the sake of completeness. But, as in [1], so here we have to refer those who 
would like to see the full calculations to the manuscript which is deposited at 
the Institute for Advanced Study. 

Finally we should like to thank Mr. Lewison for his critical study of our 


previous papers, and Mr. Schild for a careful and critical reading of this 
manuscript. 


2. Notations: the gravitational equations. Since in the greater part of our 
work, we shall have to separate space and time, our notation will not be the usual 
four-dimensional one. We make the conventions: Latin indices take the values 
1, 2, 3, and they refer to space co-ordinates only. Greek indices refer to both 


space and time, running over the values 0,1,2,3. Repetition of indices implies 
summation. 


The expression 


Oy» 


(2.1) Zwie etc. stands for —= etc. 
ax” 


At infinity the gravitational field takes the Galilean values »,,, that is: 


(2.2) o_o San; 10m = 0 > No= A. 
We write: 
(2.3) Lae = Not hy; C= + kh”, 


where h,, represents the deviation of space-time from flat space and it is not 
assumed to be small. 


The h” can be calculated as functions of h,, by means of the relation 


(2.4) Lueck’ = 5,”. 


It turns out to be convenient to replace the h’s by y's which are their linear 
combinations: 


(2.5) Yur = My — 4 40 he, 
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or more explicitly: 


(2.6) yo = 4 hoo + } hss 
(2.7) Yon = hon 
(2.8) Yan = han — 4 Smnltes + 4 5 mnltoo. 


This replacement is, of course, not very important but it does simplify the 
calculations. 

Thus we can, throughout, replace the h’s by the 7’s. The equations of the 
gravitational field for empty space, 


(2.9) R,, = 0, 

can be written (see A.2) in the following way: 

(2.10) Doo + 2Anco = 0 

(2.11) Don + 2Aon = 0 

(2.12) Onn + 2ZAmn = 0, 

where: 

(2.13) Doo = — Yoojes 

(2.14) Dom = — Yomisst Yos|sm 

(2.15) Dan = — Ymniss + Ymsinst Ynsime— SmnYrejre 
and: 

(2.16) 2Ac0 = Yarjer+ 2A’ oo 

(2.17) 2Aom = Ymsjso— Yoojmot 2A’ om 

(2.18) ZAmn = — Yomjon — Yonjom + 2dmnVoe)0s 


+ Ymnjoo— SmnVooj00+ 2A’ mn - 

In these formulae, all the linear terms are written out explicitly, while A’,, 
stands for all the non-linear terms in the y's. The division of the linear expres- 
sions into those belonging to #,, and those belonging to A,, may seem artificial 
at thismoment. In anticipation of further development, we shall remark here, 
that, in the actual approximation procedure, by which we shall solve the gravi- 
tational equations, these linear terms collected in A,, will behave like the non- 
linear terms. 


3. Lemma. We mentioned in the introduction that the differential equa- 
tions of motion will be derived by forming surface integrals. The technique 
of calculating such surface integrals will reappear many times in this paper and 
it is based on a lemma to which we shall refer asthe lemma. Here we shall give 
its formulation and its proof. 

We have a set of functions: 

(3.1) Tr. 

It is immaterial whether these functions of x* have tensorial character, or not. 
The bracketed indices are Greek, or Latin, and they will not play any role in 
our argument. But we do assume that these functions are skew-symmetric 
in the indices k, /: 

(3.2) Fy... yar — Fe... 
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We now form an integral 
(3.3) | Fi. ° ayia, dS 
(S2) 

over an arbitrary two-dimensional closed surface that does not pass through the 
singularities of the field. In (3.3) 

> 
(3.4) n= cos (x*, n) 
are the components of the “normal unit” vector to the surface. The words 
“normal,” and “‘unit’’ are used in the conventional sense to designate the 
corresponding functions of the co-ordinates, which are implied by these terms 


in Euclidean geometry. They have nothing to do with any particular metric. 
Our lemma is: 


(3.5) | Fy. . .erpmtedS = 0. 
(Ss) 


We see that the integral (3.3) is certainly independent of the shape of the 
surface, because 


(3.6) Fi. -)rle = 0, 


and because of Green’s theorem. We can also write the integral (3.3) in the 
form 


> 
(3.7) | curl, AdS, 
(S2) 


where 
Fy. --)23= Ai; Fy. -j)21= A2; Fi. --)iz = a. 
But (3.7) and therefore (3.3) can be changed, by Stokes’ theorem, into a line 
integral over the rim of the surface. If the surface is closed, the rim is of zero 
length. Therefore, our lemma as expressed by (3.5) is proved. 


4. Surface integrals. We treat particles of matter as singularities of the 
field. Let us assume p particles and the knowledge of their world lines. Thus 
we denote by 


(4.1) E(x): 5 = 1,2,3,...,P, 
the world-line of the sth singularity. Here and later, the index written on the 
top will always label the particular singularity. 

The gravitational field, that is the y's, will depend on the x*’s but also on the 
t’s and their time derivatives. The equations that the y’s fulfill are 
(4.2) $,,+ 2A,,= 0. 

At an arbitrary moment x°, let us surround the sth singularity, and it alone, 
by a closed surface. Then: 


(4.3) | “Gat 2A,n)mdS = 0, 


where the s over the integral indicates here, and later too, that the integral is 
to be taken on a two-dimensional surface surrounding the sth singularity and 
it alone. 
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We shall show that 


Indeed it follows from the definition (2.14) and (2.15) of #,, that it can be 
written in the following form: 


(4.5) Dy = Foyer 
(4.6) Fy)kt= Vatik— Vekit— Sx Virie A SyrVerir 


But F(,): is skew-symmetric in k and /. Therefore (4.4) is fulfilled. From 
it and from (4.3) we deduce: 


(4.7) | 2aamas = 0. 
Also, because of the structure of ,, we easily verify: 
(4.8) ®,ni2 = 0, 
therefore also: 
(4.9) Again = 0. 


Equation (4.9) tells us that no surface integral of the form (4.7) can depend 
on the shape of the surface. But equation (4.7) tells us more; namely, that 
such an integral vanishes. 

The 4 surface integrals in (4.7) can give us no relation between the space 
co-ordinates of the field, because the surface is entirely arbitrary. They can 
only give us relations between the co-ordinates of the singularities and their 
time derivatives. Thus we may have at most 4 differential equations. Antici- 
pating the later development, we may remark here that these equations will 
determine 3p functions of time 

£*(x%), 


that is, the motion of singularities. 


5. The method of approximation. The problem before us is to solve our 
field equations and to deduce the equations of motion. This we shall do by 
a new approximation procedure. Let us assume a function ¢(x*, A) developed 
into a power series in the parameter \ (for small \): 

@o 
(5.1) o(x",r) = No + No + Not--- = F Ao. 
0 1 2 i<0 1 

The indices below indicate the order (i in \' is always the exponent, not 
the index). 

If the function ¢ varies quickly in space, but slowly with x°®, then we are 
justified in not treating all its derivatives in the same fashion. The derivatives 
with respect to x° will be of a higher order than space derivatives. We can 
formalize the procedure by introducing an auxiliary time r, 


(5.2) r= x, 











wa ~~ 
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so that derivatives with respect to r can be treated on the same footing as the 
space derivatives: 


d¢ dg 
(5.3) =<—— = —,X = AY, 0. 
Vie ax® ar 
We conclude: the “stroke differentiation” of a quantity with respect to x®, can 
be replaced by the ‘comma differentiation” with respect to r if the power of A 
with which this quantity is associated is simultaneously raised by one. To ex- 
press this explicitly we use numbers under zeros, written after the comma, e.g.: 


(5.4) 7 nn 0 = AUF An ,0 Or: AZ en \00 = AUT On , 00+ 

21 a 1 2 2 2 
From now on, all differentiations will be with respect to (r, x', x*, x*) and they 
will be denoted by commas: 


(5.5) el el Oe 


Thus we shall develop all functions that appear in the field equations in 
power series in A. We start with the 7’s in the following way: 


Yo = Myo0+ A*yoo+ A*yoot+: - - 
4 6 
(5.6) You = Myom+ Ayom+: * 
5 
Ymn = AYymnt A*Ymnt: > 
4 6 


Why do we start with different powers of \? This is an assumption, but it 
can be justified heuristically. Assuming for a moment the usual energy 
momentum tensor for matter, we have, for a quasi-stationary field, approxi- 
mately: 





Ayoo = — 2p 
dx™ 
AvYom — = 2p ommm 
(5.7) dr 
A _ 9 dx™ dx" 2 
—_ , aa 
therefore 
(5.8) Ymn™ AYom™ d*00, 


and it is pure convention that we start with \? for yoo. 

The other question suggested by (5.6) is: why do we omit the odd powers 
of \ in the developments of yoo, Ymn, and the even powers in yon? Indeed, we 
could have introduced all powers in (5.6). A more thorough investigation 
shows that our choice (5.6) means that what we are doing here is similar to 
the procedure in electro-magnetic theory when we take not the retarded, but 
the half-retarded plus half-advanced potentials [3]. 
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All the functions that will appear later are gained from the y's by sum- 
mation, multiplication, differentiation. Thus to every component, the fol- 


lowing rule applies throughout: Any component having an {ode number of 


zero suffixes will have only ae powers of \ in its expansion. 


6. Field equations and the approximation method. We go back to the 
field equations 
(6.1) $,,+ 2A,,= 0 


into which we introduce the 7's in their power-series development. Thus the 
(00) equation in (6.1) can be written: 


(6.2) TA" (Hoo + ZAoo) = 0. 
i 2l 2l 


Now we cut up (6.2), and the other field equations, into equations for each 
approximation step. We write them down in the following form: 


(6.3a) Boo + 2Ao0 = 0 
2l—2 21—2 
(6.3b) Domt+ 2Aom = 0 
2i—1 2i-1 
(6.3c) Pant Z2ZAmn= 
2l 2l 


Let us analyse more closely the structure of (6.3). Remembering (2.13) to 
(2.15) we can write more explicitly: 


(6.4a) Doo = — Yoo, rr 
21-2 21—2 
(6.4b) Dom o- Yom, rr Yor, mr 
2i-1 2i-1 2i—1 
(6.4c) Onn= — Yun, rrt Ymr,art Ynr,mr— SmnYre.re 
2l 2) 2l 2i 2l 
and: 
(6.5a) 2Ao0= Yre.ret 2A’ oo 
2i—2 2i—2 2l—2 
(6.5b) 2ZAom = — Yoo, om Ymr,or+ 2A’ om 
2i-1 2i—2 1 2j—-2 1 2i—1 
2Amn = = Por Yon, m+ 25 mn or , Or 
-11 2i—1 1 2i-11 
(6.5c) 


+ Ymn,00— SmnVoo, 00+ 2A’ man. 
2i-2 2 2i—-2 2 21 
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Let us now assume that: 


(6.6a) Yoo --+ ‘Yoo 
2 21-4 
(6.6b) Yom ++- Yom 
3 21-3 
(6.6c) Ymn ++: Ymn 
4 21-2 


are all known. Then yoo can be found from (6.3a). Indeed Aoo contains only 
21-2 21-2 
terms already known, since mn is known and A’ oo is non-linear and can there- 
21-2 21—2 
fore depend only on the known y's. The same is true for (6.3b) and (6.3c). 
The unknown functions are contained in #’s; the known functions in the A’s. 
The oo, already found from (6.3a), appears as a known function in Aom. Simi- 
21-2 2-1 
larly yom found from (6.3b) appears as known in Am». Indeed we see now the 
-1 21 
reasons for our division of linear terms. 


Thus our equations (6.3), if solved, will give us 


(6.7) Yoo, Yom, Ymn» 
2-2 2-1 2 
and if such a procedure converges, we can determine the field to any approxi- 
mation we wish. 
The important question to consider is: are the equations (6.3) always 
solvable? 


7. The divergence condition. We go back to our equations (6.3). The 
first of them, that is 
(7.1) Poot 2Aoo = 0 
U-2 2-2 
is, because of (6.4a) and (6.5a), a Poisson equation, where Aoo is known. There 
is no difficulty in integrating this equation and finding yoo. Next we have 


21-2 
(6.3b), and because of (6.4b), we see: 
(7.2) Dom, m= 0. 
21—1 
Thus the next three equations can be integrated only if 
(7.3) Aom,m= 0. 
21—1 
But Aom is already known. Therefore we must be sure that our procedure 
2I—1 
leads us to Aom satisfying (7.3). Similarly the last six equations (6.3c) lead us 
2—1 
because of 
(7.4) Pan, n= 0 
21 


to the integrability condition: 
(7.5) Amn, n= 0. 
2 
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We shall prove that (7.3) and (7.5) are satisfied, if the field equations are 
satisfied in all the previous approximations. 


The tensor 
(7.6) Gy = R,, — $2,R 
satisfies the Bianchi identity 
a 6 _ B — 
(7.7) Gin + fk G, {sy Gg 0. 


We assume that all field equations up to the order (2/ — 2) are satisfied, 
that is including 


Poo + ZAoo = O. 
21-2 2i—2 


We know, that putting $,,+ 2A,,= 0 is equivalent to putting R,,= 0. From 
A.2 follows: =. 

(7.8) $+ 2A= — 2(Ryo— $ 100 Rag), 

which means, that our $,,+ 2A,, are a linear combination of the R,,. Thus, 
if our field equations are satisfied, then we have: 


(* = Go = --:- = Go = 0 
2 4 2i-—2 
~ G nn = G nn = : = G nn = 0 
(7.9) 3. 5. 1-3 
Gan = Ginn ™-*:- = Gmn = (0. 
2 4 2i-—2 


Let us write down the zero Bianchi identity of the order (2/ — 1). From the 
left-hand side of (7.7) we have, putting » = 0, the following linear terms: 
(7.10) — Gom, m+ Goo,o- 

21-1 a—2 1 
The non-linear part contains the products of the G’s and the 7's. But because 
of (7.9), both the non-linear part of the Bianchi identity and the second expres- 
sion in (7.10) vanish. Thus the zero Bianchi identity, together with the field 
equations give: 


(7.11) Gom, m= 0. 
2-1 

Because of (7.8), (7.6) and (7.2) this means: 

(7.12) Aom, m= 0. 
2i—1 


Going on to the next approximation step, let us now assume that besides 
(7.9), we have also: 
(7.13) Gom= 0. 
2—1 
Putting into Bianchi identity (7.7) » = m, we have in the 2/ order, because of 
(7.9) and (7.13): 
(7.14) Gan, n= 0 
21 
and therefore because of (7.4), (7.8): 


(7.15) Ama.n= 0. 
2i 
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Thus the divergence conditions are satisfied in each approximation step, 
though not identically. They are satisfied because of the Bianchi identities 
and because of the previous field equations. 


8. The surface condition and the equations of motion. We now approach 
the most essential part of our argument. We are faced with the task of solving 
the following system of equations: 


(8.1a) Poo ~_ 2Aoo = 0 
21-2 2—2 

(8.1b) Dom a 2Aom = 0 
2i—1 2-1 

(8.1c) Gan + 2ZAnn = 0. 
21 21 


We know that because of the Bianchi identities and because (as we assumed) 
similar equations had been solved in the previous approximations, we have 


(8.2) Aom, m= 0; Amn, a= 0. 
2i—1 21 


Let us also remember, that there is no difficulty in solving (8.la) which is a 
Poisson equation. But what about (8.1b) and (8.1c)? 

Before we return to this fundamental question, we wish to discuss the start 
of our approximation procedure which determines the character of our cal- 
culations. 

In (8.1) we put / = 2 and write the first two equations explicitly: 


(8.3a) | a ae = 0 


(8.3b) — Yom, ss + Yos,ms= Y00, Om- 
3 3 2 1 


The character of the entire solution will depend on the choice of the harmonic 
function we take as the solution of (8.3a). As we are interested in solutions 
representing particles, we shall write: 


ss 
ie 2¢;¢= bs {— amy} 
(8.4) \s s s s 
lv = [(xt— g*)(x* — gh) P= (7). 
Here ¢ is the “distance” in space of a point from the sth singularity. 
We leave it undecided, for the moment, whether m is a function of time, 


or aconstant. Now we introduce this yoo into (8.3b) and again obtain three 
2 
equations for the three functions yom. But is (8.3b) always solvable? True, 
3 


the divergence of both sides vanishes. But this is not sufficient. The surface 
integral of the left-hand side of (8.3b) vanishes, as follows from the lemma. 
But then the surface integral of the right-hand side of (8.3b) must vanish too. 
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If we calculate the surface integral around each singularity, we find (see A.4) 
that it vanishes only if 


d ; s s 
(8.5) a (m) -9.0.°F= 0, 


that is if the m’s do not depend on time. This is so, because 


s s s e k 
(8.6) Vio= —v nk; (= #) 


and because only expressions proportional to r~* can give a contribution to 
the surface integral. ‘Thus, going back to (8.4), we have to assume that 


2 
(8.7) m,m,m,...,m™ 
ss & 2 
are constant. 
These constants (8.7) can be positive or negative. We shall assume that m 
are positive. Indeed, by taking the first particle and removing all others, we 
see that m is its gravitational mass, since for large r the field is that of a particle 


with gravitational mass ™ This is the same constant of integration that 


appears in the Schwarzschild solution, since our field for one particle is that 
of a Schwarzschild singularity when r is large. Thus we shall have to exclude 
from our solution negative gravitational masses. But then we must also exclude 
dipoles and poles of higher order. 

Yet if we try to solve (8.1) we see (the details will be presented later) that 
we cannot do so without adding certain poles and dipoles to Yoo. This we shall 


have to do, in order to insure the integrability of (8.1) in each approximation. 
But then the solution of the total field will contain dipoles which are not 
allowed, since they represent physically meaningless solutions. We shall have 
to remove them after the total field has been calculated. This can be done by 
restricting the motion of particles. That is, the condition that the dipole field 
vanishes will give us 3p ordinary differential equations for the motion of p 
particles. Thus the motion is undetermined in the approximation procedure. 
It becomes determined after the approximation procedure is finished and the 
dipole fields are removed. 

In practice, we find solutions both for the field and for the equations of 
motion only to a certain approximation, say 2n. We obtain the equations of 
motion to the 2m approximation, by removing all the dipole fields to such an 
approximation. 

Although we have developed our field equations with respect to an arbitrary 
parameter \, this \ can be absorbed by the actual equations of motion through 


the change of scale in m and 7, so that \ is absent from the final form of the 


equations. 
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We have given a general outline of our treatment. Turning to the details, 
let us see why (8.1) will not, generally, be integrable. We know, from the 
contents of Sec. 4, particularly from (4.4) that the surface integrals of the @ 
functions vanish. Although this was proved for the total field it is equally 
true in each approximation step, since the proof made use only of the structure 
of the @’s, which is the same for the total field, as for the field in each approxi- 
mation. Thus we have: 


Ss s 
(8.8) | $on,dS = 0; | $,,,n,dS = 0. 
21-1 21 


But then our equations (8.1) can be self-consistent, only if we have: 


(8.9) | 2Aon-dS =0; | 2Amn,-dS = 0. 
21-1 21 


But the A’s in (8.1) are already known; they are functions of the known field 
calculated in the previous approximation steps. Therefore we can calculate 
the integrals (8.9) and find whether they vanish or not. 

At this point it is convenient to introduce a new notation. Because of (8.2), 
the surface integrals (8.9) will not depend on the shape of the surface, but only 
on the singularities and their motion. Thus the surface integrals, even if they 
do not vanish, can be functions of r only. 


We write: 
(8.1 ; : ; 
(8.10) 2 | Qhen dS = Colt) = Co 
ym 1-1 21—1 2-1 
(8.11) 1 


s Ss Ss 
— | 2AmentrdS = C,(r) = Cy 
4r ai 2 2—1 
and assume that we have calculated the C’s. If they vanish identically, and 
if they vanish always as we proceed with our approximation, then our equations 
are self-consistent. 

Let us assume, however, that the C’s in (8.10) and (8.11) are mot zero. Then 
(8.1b, c) cannot be solved. There is nodifficulty in solving (8.la). This equation 
is of the form 
(8.12) Yoo,rr= 2Aoo, 

a—2 21-2 
where the right-hand side is known. We see that the solution of this equation 
is determined only up to an additive harmonic function. Thus we can add to 
any solution either single ‘‘poles” or “‘poles’’ and “‘dipoles.”’ 

By adding single poles we can insure the integrability of (8.1b). Then by 
adding dipoles we can insure the integrability of (8.1c). We could have done 
all that in one step, adding poles and dipoles, but the division into two steps 
makes for a simpler presentation. 


s Ss 
After finding yoo from (8.12), we calculate Cy» and, in general, find Co# 0. 
21-2 u—1 21 
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We then replace in (8.1b): 


(8.13) yo by vo — LD 4my 
2-2 2-2 s 21—2 


$s 
where m are certain functions of time to be determined soon, and y’s are the 
21—2 


functions defined in (8.4). Of course this change in yoo induces a change in 
21 


s -2 
Co. Indeed, 
2—1 
2Aom changes now to 
(8.14) wie 


|2Aon-+ “(4 m v), Om) 


1 s 2i—2 1 


as follows from (6.5b) because yoo appears in Aom only as — Yoo,mo- Now 
21-2 2-1 u-2 1 
obviously the old surface integral 


om RS | 2Aon-dS = Co 
4a 2i—1 2-1 
changes into A.4 
(8.16) Co — 4m, 
2—1 211 


therefore it can be made zero by choosing 
(8.17) 4m = Cy. 


Thus by adding a pole we can insure the integrability of (8.1b). The next step 
is to insure the integrability of (8.1c). Thus we assume that ‘oo, Yon are known, 


that (8.1b) is integrable and we have once more to return to oo looking for a 
21-2 


different solution of (8.la) so as to insure the integrability of (8.1c) without 
destroying the integrability of (8.1b). 
We replace now our yoo (containing the additional poles) by 
21-2 


D s Ss 
(8.18) vo — > S- vir. 
2-2 = s=1 21-2 
These are additional dipole solutions, and we assume that no other dipole 
expressions are contained in yoo. Again the S, are functions of r only, to be 
21-2 


determined later. The yoo now contain the single pole solutions so as to 
21-2 


enforce the integrability of (8.1b). We can easily see what change in yom is 


induced by (8.18). The answer is, that yo, changes into 
21-1 


(&.19) Yom— DL (Sm ¥),0- 
2l— 1 


1 s 21-2 








(§ 
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Indeed, if the old yom satisfies the original equation (8.1b): 
21—1 


(8.20) Yom, ss YOs,mse= Yms, 0s 700, mo+ 2A’ om 
2i—1 2i-—1 2i-2 1 2i-2 1 2i—1 


then Yoo, Yom with the additional expressions written out in (8.18) and (8.19) 


2-2 2I—1 
satisfy the equation too. This isso, because 2A’, being non-linear can contain 
21 
neither yoo nor Yom. Therefore the addition of dipoles does not affect the in- 
—2 21— 


tegrability of (8.1b). 
Now the last and decisive step: we replace in (8.1c) yoo, Yom, by the new 
2-2 2i—1 
expressions according to (8.18) and (8.19) and adjust the S’s so that the surface 


integrals will vanish identically. This requires a somewhat more lengthy 
calculation. 


Written out explicitly, equation (8.1c) is: 
Ymn, ss zane. ns Yne,met Smn'Yre,re 
2 


so- ae on a. ont 25 mn or, ort Ymn, 00 bmn 00, oo+ 2A’ mn 
1 


2-11 2=2'2 2-23 2 
(8.21) 
= 2Aun- 
21 
We introduce into (8.21) 
(8.22) yo - LD S++ 
U=-2 3 2-3 
(8.23) Yom— LD (S. ), 0 
g—-1 s \ai-2 / 1 


for the old yoo, Yom. We now obtain new expressions added to the old 
Amn. The difficulty is, that now the contributions come not only from the 
21 

linear expressions, but also from A’m, which will contain terms of the type 


Yoo - Yoo. The result of the calculations is given in A.8, and contains many 
—2 


amnatuen of which we shall here write only the first three which arise from 
the linear terms (the others, as we shall see, are unimportant). Instead of the 
old 2Am, we have: 

2 


ZAmn 
21 


(8.24) + y (5. at Sut ea bmn d SH, .) 
+... 


where the dots at the end indicate the omitted expressions. As we are here 
discussing the problem of surface integrals, we are justified in omitting them 
because they do not give any contribution to the surface integrals. We see 
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too, that the expressions written out here have a vanishing divergence, and 
this is true for the omitted termsalso. Calculating the surface integrals (A.4), 
we find that the old surface integral 


(8.25) 1 | . _ 


changes into 


Ss s 
(8.26) Cu— Sm. 
21 21 
Therefore it can be made zero, by choosing 
s 
(8.27) Snu= Cn. 
21 21 


Thus we can always, by adding dipole solutions in oo, force the surface 


integrals to vanish identically. 

By proceeding in this way, we accumulate single poles and dipoles, and the 
additional expressions in oo are: 

p $ s s s 
(8.28) — ENE (tm +S, v,). 
i s=1\ 2I-2 21-2 

We violated our rule of not introducing dipoles. However, this was done 
for yoo only. Wecan, at the end of the approximation procedure, annihilate 
all these additional dipole expressions by taking 


s 
(8.29) > a**s, = 0. 
i 21-2 
Differentiating this twice, we obtain, because of (8.27): 
Ss s 
(8.30) Er A*S,.= ¥ A"Cr= 0. 
7 21 i 21 


These are the 3p equations of motion. Thus the motion is determined, if 
dipole solutions are rejected. 
On the other hand, the m’s can be calculated from the Co's according to 


(8.17). Denoting the total coefficient at y by — 4M, we have: 


s Ss s s 
(8.31) M = Xm + Mm + Nm +: -- 
2 4 6 
Ss Ss Ss 
where m, m,. . . are functions of the original constants m and of known func- 
4 6 2 


tions of the time. 
The equations (8.30) and (8.31) will contain only a finite number of terms 
depending on the order to which we wish to carry out the actual calculations. 


9. On the choice of a co-ordinate system. We shall now see that it is 
possible to simplify our equations through the proper choice of a co-ordinate 
system. Let us assume that 


(9.1) 1* 00, Y¥*om > ¥*mn 
2l a 





ar 
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are solutions of our system (6.3), where the #’s and A's are defined by (6.4) and 
(6.5). Then we can show that any 


Yo = Y 00 


2i-—2 2i-—2 
(9.2) Yom = Y*om + Go, m 
21-1 21-1 21-1 


Yan = Y*mn + Om,n + Qn,m — Smn Or, rt bmn Qo, 0 
2i 2 21 24 2! 2-11 


with ao, a» arbitrary are also solutions of our equations. This can be shown 
—1 2 


just by straightforward substitution in (6.4). A simple calculation shows that 
all the a’s vanish from these equations. Thus we can, at each approximation 
step, impose four conditions upon the field. Let us choose, as is usually done, 
the following four co-ordinate conditions: 


(9.3a) Y00,0— Yor,r= 0 
21-2 1 2i—1 

(9.3b) Yom,0— Yar,r= 0. 
2i-11 2 


Indeed, if y* do not satisfy such a condition, then a’s can be found that 
ensure it. The equations for the a’s are: 


(9.4a) ao ,rr= ¥*00,0— vor. 
2i—1 2i—-2 1 2i—1 

(9.4b) an ,rr= 1*om,0— Yar, 
21 2—1 1 2l 


With the co-ordinate condition (9.3) our system of equations is considerably 
simplified. Equations (6.3) now become: 


(9.5a) Yoo, rr = Yoo, 00 + 2A’ oo 
2i—2 2i—4 2 2i—2 

(9.5b) Yon rr = Yom, 00+ 2A’ om 
2-1 2-3" 2 21 

(9.5c) Ymn ,rr= Ymn, oot 2A’ ans 
2l 2-2 2 2 


which together with the co-ordinate conditions 


(9.6a) 00,0 — Yor,r = 0 
2i-2 1 2i-1 

(9.6b) Yom,0—~ Yar,r= 0, 
2i—1 1 2l 


now form a symmetrical system of equations, where in (9.5) all the known 
functions on the right-hand side are at least two orders lower than those on 
the left. 














226 A. EINSTEIN AND L. INFELD 


The surface integrals that must vanish and which give the equations of 
motion are: 


(9.7a) | Yom,00— ‘Yoo, om+ 2A" \’on) tmndS = 0 


2I1—3' 2 21-21 


s 
(9.7b) | Ynm,00— Yno, om + 2A’ am) Nm dS = 0. 
2—-2°2 2-11 

We can deduce them from our old formulae, using the lemma, or directly, 
differentiating (9.6), adding to (9.5) and using the lemma. 

If, as in Sec. 8, we now introduce dipoles in order to satisfy (9.7b), we do 
not violate (9.6a). 


Sometimes it is more convenient to use other co-ordinate conditions. For 
example, the one used in the actual calculations is: 


(9.8a) Yoo,0— Yos,..= 0 
2I—2 1 2i-1 


(9.8b) Yma,n = 0. 
21 


The equations then are: 


(9.9a) YO,rr = 2A’ oo 
21—2 2i—2 
(9.9b) Yom,rr = 2A’ om 
2i—1 2i—1 
(9.9c) Ymarr = — Yom,0n— Yon, ont bmn Yoo , 00+ Ymn, 00+ a mn 
2l 2i-1 1 2i—1 1 2i-2 2 2i—2 2 
_= ZAmn 
21 


and the surface conditions are: 


(9.10a) | 2A’ om — 700, om mds = 0 


2-1 89. i=2 1 
(9.10b) | ZAnmtmdS = 0. 
21 


The question arises: to what extent does the co-ordinate condition influence 
the equations of motion? We shall return to this problem in the last section 
and we shall show that the equations of motion to the sixth order do not depend 
on the choice of the co-ordinate system. 








( 


of 


or 
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10. The Newtonian approximation. We shall discuss now the first three 
equations for / = 2. The equations are: 


(10.1) Yoorr = 0 
(10.2) Yom,rr = 9 
3 
(10.3) Yanm,rr = ZAnm- 
4 4 
The co-ordinate conditions that we accept are: 
(10.4) Yor,, —~ Yoo = 0. 
3 , 7 
(10.5) Yar, = 0. 
4 


The explicit form of Ama» is given in A.10. 
4 


The character of our entire solution will depend essentially upon the choice 
of the harmonic function we take as the solution of (10.1). As we are interested 
in solutions representing particles, we shall write: 


? ss 
Yoo= 29; ¢=2> { _ 2mi} 
2 


s=l 


i-[e-e-H]"- 


From (10.2) we see that yom is a harmonic function too, which must, however, 
3 


(10.6) 


satisfy the co-ordinate condition also. From (10.4) we have: 


5 s 8 
(10.7) Yor,,= Yoo= — po {4m dt . 
3 2 1 s 2 1 


s ° ° . ° ° ° 
The constant m, which we identify with the gravitational mass of the par- 
2 


ticles is assumed to be positive. Therefore the exclusion of dipoles, together 
with the field equations and the co-ordinate condition determine uniquely Yon: 
3 


sss 
(10.8) Yon = L 4myi". 
3 s 2 1 
To this yon we could add, according to (9.2) the gradient of any function 
3 
and in this way obtain a general solution. But as our entire procedure consists 


s 
in employing only rational functions of = er), any such addition would 
introduce new singularities (not of the character of a single pole), or a non- 
Galilean field at infinity. Thus we should regard yo, in (10.8) as character- 
3 


izing the problem of particles, regardless of whether we introduce the co- 
ordinate condition (10.4) or not. 
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Just for the sake of simplicity, let us now restrict our consideration to two 
particles and write (omitting the indices below m, ¢, f, 2g): 


ge=f+g 


11 22 
(10.9) f = — 2my ; g = — 2my 





1 2 
L i = n’; t= of 

The next step then, since the surface integral (9.10a) vanishes for / = 3, 
because 


(10.10) | (24’en — yo19n) NmdS = -| 100,0m %m dS = 0, 
3 2 21 


1 
is to determine 


1 if} 
Ca z| 2A mr mr dS 
4 


(10.11) aj 4 
2 1 2 
C= if 24.245. 
4 4a 4 


If we wish to finish our approximation procedure here, the equations of 
motion up to the fourth, or as we shall call it, the Newtonian approximation, 
are: 


1 2 
(10.12) Cn= 0; C= 0. 
4 4 
All we have to do now is to calculate the surface integrals, according to the 


method outlined in A.4. The result of this particular calculation is given in 
A.10. It is: 


| 
—) 


: 3 (. o. 
(10.13) 4 Cm(r) = sm {em of - - 


Zim = £,m for x*= 7° 





fim = fim for x*= £*. 


The form (10.13) is actually independent of the variables x*. In the last 
equations we see that Z,m, say, is obtained by differentiating g with respect 
to x* and then by replacing x* by 7*. But the result will be the same if we 
first replace x* by n° and then differentiate with respect to n* or ¢*. Thus: 


_ 8) _ _ aetr) 


dn* o¢* 


22 
(10.14) g(r)= — = r= (n*—f")(n*— f°). 





a+ s+ —_— © @® 
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We can, therefore, think of our equations of motion as involving the differ- 
entiation of functions depending only on the position of singularities, as is 
characteristic of the theories based on the concept of action at a distance. 
Indeed, we see that our equations are precisely the Newtonian equations of 
motion, deduced here as the first approximation from the field equations. The 
treatment of p particles (instead of two) does not add any new difficulties if 
we deal with the Newtonian approximation only. 


11. Transition to the next approximation. We wish to go now beyond 
the Newtonian approximation. But then we must calculate yma, since Amn 
4 


6 
depends on mn. The characteristic feature of this method is that generally, 
4 


if we wish to find the equations of motion to the 2/ approximation (inclusive) 
then we do not need to calculate Youns because Cm does not contain it. But 


now, if we wish to go one step darthine we must find Yonn for which the equations 
are: 
(11.1) Ymarr = 2Amn- 
4 4 

This is “the transition step” that we have to take before proceeding to the 
next approximation. These equations are integrable only if we do assume 
Newtonian motion. Otherwise we would have to add dipoles. Yet if we wish 
to proceed only to the next approximation we may assume Newtonian motion 
and additional expressions induced by the dipole fields are not necessary. 

If in (11.1) we assume Newtonian motion, then (11.1) can be integrated, 
because the surface integral of A,» vanishes then. But if we do this, we intro- 

4 


duce Newtonian motion into Am». This is admissible because any difference 
6 

between A calculated this way and A calculated with the proper motion is of 

order A. * Thus since we do not sicaitiie to go beyond A we may ignore the 


additional dipole fields. It is for this reason that the nnehens special calcu- 
lations in [1] were correct, but the general theory was not. 
We shall now solve 


(11.2) Yma,rr = 2Ann = — Yom,0n —~ “Yon,0m + 25 mn ¥,00 
4 4 3 1 3 1 


-_ 2¢¢,mn — OmP,n + 3 bmn a 


assuming the Newtonian equation of motion, i.e. (10.13). 

We can ignore the dipole expressions because we are interested only in the 
equations of motion to the next approximation. But, for the same reason, we 
are interested only in those expressions in ym», which give a contribution to the 

4 


corresponding surface integral of Amn. 
6 
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An inspection of Am» (A.12) shows that we need only the knowledge of yma 
6 4 


in the neighbourhood of the singularities, and we may ignore in it the terms 
which do not become infinite as 7 + 0, since the surface integral due to these 
terms must vanish (see A.12). On the other hand y,, which also appears in A 

4 6 


should and will be calculated in the entire space. 

In the equation (11.2) we have, on the right-hand side “cross products,” that 
is, products belonging to different singularities. Because of them (11.2) can 
only be integrated in the neighbourhood of the first singularity, say. The expres- 
sion arising from the second singularity can be expanded into a power series 
near the first singularity. Retaining all the expressions that may give some 
contribution to the surface integral and those only, we have in the neighbourhood 
of the first singularity: 
[Ymn = { f[(x*— a) i" + (x™— 2) 9" — Smn(x*— n*)i'}}, 0 
+ { gl(x*— E+ (e™— EI" — Snnlx*— FET, 0 
(11.3) } + d Pf mfin + 7 7g m Bn 
— fm(x"— n")g 
\ + Omnf + Binng.- 

Here only the expression 

- S.m(x*— n")Z ; g = for x*= n*, 

is due to the interaction terms. The two last expressions are the additive 
harmonic functions (dipoles are excluded) and they are determined by the 


co-ordinate condition 
(11.4) Ymr,r = 0. 
4 


The result is: 





(11.5) Bmn= 2™{"+ Smnf, 
F=f): @ = g(r); P=" — Fn" f°). 
But, let us say once more, that all this is true only if the Newtonian motion 
is assumed. 
Finally, as we mentioned before, y,, can be calculated rigorously. The 
result is: : 


(11.6) a 2m 7.00 — 2m 700+ i¢ 
+ af + Bg. 


Here the a and @ are determined so that near the singularity (11.6) will be 
consistent with (11.3) form =n =r. The result is: 


(11.7) a = 27°7'+32 
B = 2+ of. 


Thus our transition steps are accomplished. 


= = 259" + 5 mnZ 











‘ve ‘" 
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12. Beyond the Newtonian approximation. We write down the next field 
equations: 


(12.1a) yor ™ hn = 3 Yr er 
(12.1b) Tomer = 2. om = Ys Y0s,m— Yam Yes — 3¢,00,m 
(12.1c) Ymayrr = ZAmn- 

6 6 


The explicit expressions for A,,, are quoted in A.12. The solution of (12.1a) 
6 
is simple: 
11 22 
(12.2) yo = — £¢’— 4myp — 4m. 
4 4 4 

As we know from the general theory, the arbitrary harmonic functions have 

to be determined in such a way as to make (12.1b) self-consistent, that is, the 


corresponding surface integral must vanish. 
The co-ordinate conditions, are here, as before, 


(12.3a) Yor,r —~ Yooo= 0 
5 4 1 
(12.3b) Ymr,r = 0. 
6 


Because of this, the conditions for solvability of (12.1b, c) are: 


(12.4a) i | }2A'om— 700, om} tim dS = 0 
4n 5 4 ‘1 

(12.4b) = | 2hermeds = 0. 
dn 6 


We have in (12.4a) the equations that determine m. The result of evaluating 
4 


the surface integrals in (12.4a), (see A.12) is: 


r 


ma sid atirt 3 \ 2 (shirae— shin + ) 

ek ak See od let | weit 
. os Relies 17h 2 (aise i * ) 

(12.5) | natin} e+ os “% me*t mm ~ 
m= m; m= m; t= (9*— f*)(n*— £°). 





The next step, after the self-consistency of (12.1b) has been insured is to 
calculate the yo,. We need them, because they enter into the next surface 
5 
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integral. Including only relevant terms that can influence the surface integral 
we find near the first singularity: 


yom= — Ei finhe T+ ELT 
+ § (x*— 0°)(i*— £*) fem 
(12.6) — (x™— 9™)fe..(i°— §*) 
+ } (x*— a )fmo* {2 + Z,-(x"— 9)} 
+ & (x*— 0°){ fo." + fbi } 
+ aomf. 





Again aom is determined from the co-ordinate condition (12.3a) and the 
result is: 


(12.7) @om= — 9° 9° G+ Bi" — Bi”. 
Now the scene is set for the last and most difficult calculation: 


1 1 1 
(12.8) C.= — | ZAmntmdS. 
6 6 


Some remarks about this calculation are made in A.12, and partial results 
given. We obtain: 


1 , ; 2 1 1 
Ca= — ad. 2. {| var 3 rere 45%? —4— = | aS (-) 
6 2 r r | dn™ 


-s(-m_ om “ms Aveem 2 1 1 = err 
+ [at(im— om) + 3ari+— afi] (7)+; tit. 


Thus the equation of motion belonging to this stage of approximation is: 
1 1 
(12.10) A‘C,, + ASC. = 0. 
4 6 


. _— : 1 2 
We can now re-absorb the \’s by substituting new units for r and m, m: 
old r = \- new 7; old mass = \~*- new mass. 


Preserving the old symbols for the new units we have for the equations of 
motion of the first particle: 


\ . 2 1 
i"—m a = in {| iran § it — 4qtje— 4 - * | * (1/7) 
” r r jon 
: ‘ — 0 
(12.11) + [4arGm— a) + Bari 4b] (/) 


+ 1 > cer \ 
2 dn*dn’dn™ J : 
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The equations of motion for the other particle are obtained by replacing 


m, m, n, f, by m, m, [,”, 
respectively. 

These are the equations of motion of two particles. They can be integrated 
and conclusions concerning perihelion motion of a double star can be drawn 
from them [5]. The entire method can also be adapted for the case of a 
charged particle in an electromagnetic field [4]. 


13. The equations of motion and the co-ordinate condition. The con- 
tents of the last three sections are not new. Its presentation, however, is 
different than that given before in [1] and [2], since it has been adjusted to the 
new theory. There is one more question that we wish to answer and which 
we did not treat before. It is possible to do so only now after the general 
theory has been perfected. We ask: To what extent do the equations of 
motion as formulated in (12.11) depend on the particular choice of the co- 
ordinate system? 

We reject any particular choice of co-ordinate system and write the first 
two equations: 

(13.1) Poot hn _ ‘Yoo,er = 0 


(13.2) Dom t+ 2ZAom= — Yom,rr+ Yor,mr— Yoo,mo = 0. 
3 3 3 3 a Ss 
We assume that we start our approximation procedure with the same ‘oo 
2 
and ‘yom functions as we did before. But from now on, while dealing with the 
3 


rest of the equations we shall look for general solutions not restricted by any 
additional co-ordin: te conditions. 
Thus the equativ: 2 that we wish to consider now are: 


(13.a) nat ZAmn= 0 
4 4 

(13.3b) Poot 2Aoo = 0 
4 4 

(13.3c) Domt 2Aom = 0. 
5 5 


In the previous three sections we solved these equations, using special co- 
ordinate conditions. Let us mow call the special solutions that we obtained 
there: 


(13.4) 1" as ’ * oo, 1*om- 
4 4 5 


Knowing them, as we do, we can find the general solution of (13.3). The 
procedure is similar to that outlined in Sec. 9, only slightly different, because 
we have now a set of equations of order (2/), (21), and (2/ + 1), whereas before 
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we had a set of order (2/ — 2), (2/ — 1), and (27). Buta straightforward sub- 
stitution shows, that because of the linear expressions in (13.3), (and they alone 
enter the argument), the general solution of (13.3) is: 


(13.5a) Y¥mn = ¥* mat Guat anim Sun @r,r 
4 4 4 4 4 
(13.5b) Yoo = Y*oot a,,, 
4 + 
(13.5c) Yom = Y*om+ do,mt+ Am, 0 
5 5 5 4 1 


where a, are arbitrary. The question then is: If we substitute these new 
expressions into the A’s do we change the integrals 


(13.6) | Amen, dS, | Aon, dS, | Amn, dS ? 
5 6 


4 


As far as the first two integrals are concerned the answer is easy; A is not 
4 


changed; only linear expressions in A are affected, but the surface integral of 
5 


the additional expressions disappears because of the lemma. But it is different 
with the third surface integral. In A new terms appear containing the a's. 
6 


They appear both through the linear and the non-linear expressions. But 

these additional expressions—quoted in the last appendix—are such that their 

surface integral vanishes. Thus in the sense explained here the equations of 

motion do not depend on the choice of the co-ordinate system. This depen- 

dence would appear probably in the next approximation steps (A), but it does 
8 


not enter into the surface integral of A. This is a satisfying result, because it 
6 


is difficult to see the meaning of our co-ordinate conditions 


Yar,r= 0 
4 

(13.7) Tér,.e— Foo,e* 0 
5 4 1 
Yar,r= 
6 


and it is good to know that our equations of motion are independent of it. 
This result is general. If we have a system 


Onn + 2Amn = 0 
21 21 


(13.8) Doo _ 2Aoo = 0 
2l 2l 
Dom ~ 2Aom = 0, 
2i+1 21+1 


then the surface integral of Am» is independent of the co-ordinate conditions in- 
21+2 


troduced in this particular approximation stage. This is so, because the a’s 
21 
combine with the ¢’s in the same way in each approximation step. 








— 
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APPENDICES 
A.2 
The field equations are: 


a0) Rm NOH Hohe ted od ~ lot tok 


Introducing here the h’s as defined in (2.3) and splitting (A.2, 1) into linear 
and non-linear terms we have 


(A.2, 2a) Ro= = ; hoojsst hoswos— ; - 00 — L’ w 
Ron= = 4 honjest ; hos ne 


(A.2, 2b) + itu Shea? Vo 
Ran= —  hnajes +} hnsns + 9 basins — 3 hesins 
(A.2, 2c) + 4 hmnjoo — 4 hmoino — 4 Anoimo 


+ } hoojmn+ L' mn. 
Here L’,, are the non-linear expressions. We form now: 
(A.2, 3) — 2ARw— 4n7°Rs) = 0. 
Substituting the y's for the h’s, we see that (A.2, 3) written out is (2.10)— 
(2.18), where 
(A.2, 4) A’ y= L' yo — ¥ tye Lag - 


AA 


In calculating the surface integrals we need to take into account only expres- 
sions that go to infinity like r~*, because only such expressions will give finite 
contributions. Since all the field functions are finite (outside of the singu- 
larity), and since the contributions do not depend on the shape of the surface, 
we may ignore all other expressions. But we have to keep the surface fixed, 
because in our calculations a complicated expression whose surface integral 
does not depend on the shape of the surface, is split into partial expressions 
with non-vanishing divergence. Thus in our calculations the surface is always 
a two-dimensional “‘sphere”’ with radius shrinking to zero. Let us assume, for 
the sake of simplicity, that the space co-ordinate of the singularity is (0, 0, 0). 
We shall first give some examples of the surface integrals formed around such 
a singularity. 


Example 1. We calculate: 
0 
¥.n.dS; par; P= x'x*, 


We have: 
0 


0 2.8 8 
vn dS = — [ —_ r? sin 6dédg = — 4x. 


d r 


Example 2. We calculate: 


0 8? 
| v.n,dS = — [ ~~ rsin ddédg = — — 5,,. 
r4 3 
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Example 3. We calculate: 
[Yom na x(r)dS. 


To find such a surface integral we expand x(r) as a power series in the 
neighbourhood of the singularity: 


x = x(0)+ x,.(0)x.+.... 
The only contribution is from the second expression, that is, we have to 
calculate: 


0 ~™x* , 
x,2(0) | Vi mntnx’*dS = — x_,(0) — r? sin 6déde 


rT 
rs 


+ 3x0) [| : r? sin 6déde 
r 
8 





us 
i X,m (0) . 


In the course of our calculations we shall have to find more complicated 
surface integrals and the following table will prove to be useful: 


Table of Surface Integrals 


1 

I. 2 f at%,@dS = — 1. 
= ¥, 

er [vemas o = He... 
4n 

— . = Team = B4,. 
An 

IV. x x7 meMndS —- ys { Q5rn ine 35m Snes — 35ers bmn} . 
7 

V. » | 26 cotadS = % bss. 
An 

VI. 2 x"Y metn dS = 0. 
4n 

VII. — | x*x*Y mritn= 0. 
4x 


VIII. ~ [sx Vanteds = z bme Sir = 3 (6mt5re+ bmr Ste). 
T 


A8 


The linear terms of (8.26) give the following contribution to 2Am,: 
21 


? e -@ y @ 
(A.8, 1) Ym {50 Vat Sa Vin— bmn S, vst 1 00- 
2i—2 2-2 2 


s=1 \ 21-2 
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The non-linear terms can be found in the following way: Inspecting the 
terms in Amn(A.12, 3) we see products of yoo and yoo or, as it is there called 2¢. 
6 4 2 


Thus, if we put there the expression in (8.18) in place of yoo and write for 
4 
brevity: 
s s 
(A.8, 2) (S¥) = DS.¥ 


Ss 
we get five new terms. Thus with the abbreviation (A.8, 2) we have in every 
approximation the following additional terms: 


{ (Sm¥), nt (Sav) m— 5mn( Spy) ? } , 00 
(A.8, 3) + (SAW), rman $ ¢,n(Sr¥), rm 
+ $ ¢,m( Sr) rn % Smn?, (SW) ret ¢,mn( Sr), + : 


Only three of the linear terms give us a contribution to the surface integral. 
It is more difficult to see that the non-linear terms do not give any contribution, 
since it requires some knowledge of how to deal with surface integrals which 
is outlined in A.4, and which we shall here assume. We can write the non- 
linear terms in (A.8, 3), in the following way: 

{ , ma(Srp) } Pea { ¢, mr( Sav) } oa 
+ ¢,mn(Sn¥) + 
(A.8, 4) +  {en(S),m}.r— & e,e(Sa¥), mj. 
+ 4 ¢,r( Sa), mr 
= % Smn?, (SW), or- 

These are the non-linear expressions, and their divergence vanishes because 
y is a harmonic function. The expressions written out in pairs in (A.8, 4) do 
not give any contribution to the surface integrals, because of our lemma in 
Sec. 3. Thus the only contribution could come from the terms: 


(A.8, 5) % ¢,05nV,me— % Smn GY, 25r¥ re . 
Here only the “‘cross products’ could give contributions and we find with 
thefhelp of{the table in A.4, that the result is zero. 


A.10 
In the / = 2 approximation we have: 
[ ye= 2e = 2f + 2g 
Yon = — 2fG"— 2gt" = hon 
3 3 
hoo =g=f+g 
(A.10, 1) } h® = —ho= — ¢ 
2 
h™” = hon= Yon 
3 3 3 
han= — h™*= bmn - 
L 2 2 
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A straightforward calculation gives: 


2Aoo = 0 
2 
ZAom= — Y00,mo 
(A.10, 2) 3 : 1 
ZAmn= — Yom,on— Yon,omt 2dmn ¥,00 
4 Se 2 


jo 20¢, mn — O,mP,at 3 bmn %,29,8- 
The contributions to the surface integrals are (for the first singularity): 


— Yom, 0n => 4 mi™.4e 
— Yon, 0m — 4 mi™.4e 
mn > — $mi™.4e 
— 29¢, mn _ 3 me mAw 
— O,m%,n + —§ mE.mAm 
ZimnGs?e > 2MB.m4e 
(m = m). 
2 
A.12 
A straightforward calculation of A, A, A gives 
(A.12, 1) 2A00 = — $ 2%, ote" 


(A.12, 2) 2A’ om ¥,s Y0s,m — Y,sm Y0s — 3¢,0¢,m 
5 3 3 
ZAunn = — Yom, 0n— Yon, om+ dmnY00, 00+ Ymn,00— YY00, mn 
6 5 5 4 4 4 


— OYss,mn— $,mnYoo— P,mnVest O,msYne 
4 4 4 4 


Pins Yms— San, ar'Y or — 2¢, sY¥mn, et Y,sYme,n 
4 4 4 4 


+ + 


ined temealll 4 Y mV ssn $ 9 nVss,m— 4 Pn Y00,m 
4 4 4 


4 ¢,mY00,n + 3 bmn? .s¥rr.et 3 bmn ,s'V00,s 
4 4 4 


Vos Yon,ms— YosYom. net 27 0s'V0e,mn 
3 3 3 3 3 3 
(A.12, 3) + 5 bmnVos,rVor,s— 3 SmnVos,rV0s,r7 Yos,mY0s,n 
3 3 3 3 3 3 
+ Yom,sVY0n,s—~ $,0nYOm— 90m ont 25 mn 0s 9,08 
3 3 3 3 3 
— P0YOm,n—~ OCYon,m— O,nV0m,0— Y,mVon,0 
3 3 3 3 
2e70m,ont+ 2¢70n,0m— 25 mn¥¢ ,00 
3 3 


2009, mn— $P,mP.nt ¥ dmnG?,s,s 


++ + 


, Sb mn¥ ,0¥ ,0- 
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The surface integral (12.4a) for s = 1 is, because of (12.2), and (12.1b): 
1 
de | (¢,°Yor,m— ¢ rm On — % 09 m+ 3 PP 0m + 4(m¥) om) mdS = 0. 
5 3 3 4 


The contributions of these five expressions are respectively: 


1 
4m si - 1. : 
qdij-- = 3 B.ef*°— 4mg,.7° 
2 Sm - 7) 

( ) — = 3 £05 
(3) > = 8m Bakt+ mee" 
(4) > 2mz.. 1%" 

1 
(5) > — 4m. 

4 


Therefore: 
— 4in = me .k*+ mE. i= WME. i*— mE,.i'+ mec 
1 — aa 
— m(27* 4°+ 2) 0. 


From the last equation (12.5) follows immediately. 
The last step is to calculate the surface integrals due to A. Here a skilful 
6 


use of the lemma may save the calculation of many surface integrals. Indeed, 
2Amn can be written in the following form: 
6 


ZAmn 


(9,n¥em— ¢,s¥nm) et (GY¥me.n— PYmn,e).2 
4 4 4 4 
(8me?,r'¥rn — 5mnP,r¥rs) at (Smn?,s¥rr— Sms? .nYrr) 7 
4 4 4 4 


$ (Bmneyrr.e— 5meP¥rr,n) et (OmnVoe,0— bms'Yon,0) 0 
4 5 


a 
a 
+ 3 (Yos,mYon — Yon, m70s) et (Omn? 00s — 5ms ,0YOn) ,2 
3 3 3 3 3 3 
+ (YonYom,s— YosYom,n) et $(SmnVos,r¥or— bmseVon,rVor),s 
3 3 3 3 3 3 3 3 
ce 


(8msYor,n Vor a bmn or sor) 8 
3 3 3 3 


— Yom,on+ Ymn,00+ YosYos,mn [art+ as+ as] 
5 4 3 3 

_ 3 BmnVou,r¥ou,r—(P,n'Yom) 0 [ag+ as] 
3 

- (¢,mYon) 0+ (PY 0m,n) 0+ (PYOn,m),0 [as+ a7+- as) 

3 3 
_ 3 Smn¥ cP ,o— 4 PY ss,mn [a+ a0) 
4 
+ 5 O.nVeemt 3 P,n'700,m laut a2] 
4 

+ $ $ ,mY00,n— $ SmnP 200, lais+ ax) 
4 4 

— 99,008mn— 299,m?.n [ais+ a16) 


+ a PP 20, 00mn - [a17] 
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Because of the lemma we have to find now the surface integrals of only 17 
expressions denoted successively by ai, a2,... , ai7. The result of this calcu- 
lation is summarized in the table. Only ten types of expressions (or their 
equivalents) appear in the result. The table tells us what is the contribution 
of each of the a’s to the final result. The only a that does not give a contri- 
bution is az= Ymn.oo- 

4 


A.13 


The additional expressions in Am» induced through rejection of the co- 
ordinate condition are: ° 


2 (8mn @o,ro— 5mr@o,n0),r 
+ (G,mOn,r— O,mOr,n),+ 
+ (¢,nOm,r— 9,rOm,n),r 
+ (9,n2e,m— Y,22n,m),s 
= 2(bmn ¢,22s,r— Smr ,22s,n),r 
+ 2SmnP,er.r— dmeP,nGr,r) ,r- 
They are written in such a way, that the vanishing of each line is evident, 
because of the lemma. 
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SOME PROPERTIES OF THE EIGENFUNCTIONS OF THE 
LAPLACE-OPERATOR ON RIEMANNIAN MANIFOLDS 


° 
S. MINAKSHISUNDARAM AND A. PLEIJEL 


Introduction. Let V be a connected, compact, differentiable Riemannian 
manifold. If V is not closed we denote its boundary by S. In terms of locai 


coordinates (x*),i = 1,2,... N, the line-element dr is given by’ 
dr? = gix(x', x*,... x”) dx*dx* 
where giz (x', x*,... 2x”) are the components of the metric tensor on V. We 


denote by A the Beltrami-Laplace-Operator 
s ¢ , OU 
46 6 oe oe g a. 
Vg ax" (vie =) 


and we consider on V the differential equation 


(1) Au + Au = 0. 
If V is closed this equation will in general have an infinite number of eigen- 
values A = Am, m = 1, 2,..., and corresponding eigenfunctions ¢,,(P) where 


P isa point in V. When V has a boundary we have to consider in addition 
to (1), certain boundary-conditions in order to define eigenvalue-problems. 
We consider the following conditions, either 


(2) u=0Oon 5S, 

or 

(3) on = Q0on S, 
on 


where — denotes a differentiation in the direction of the normal of S. Eigen- 


values and eigenfunctions shall be denoted in the same way as was indicated 
in the case of a closed manifold. We assume in all cases the eigenvalues to 
have been arranged in non-decreasing order of magnitude and the eigen- 
functions to form a complete orthonormal set 


I, oi(P)de(P)dV = dix 


(dV = +/gdx'dx*... dx”). In the problem with a closed manifold and in (1), 
(3) the value Ao= 0 is a simple eigenvalue (A;> 0) with the corresponding 


1 
eigenfunction equal to the constant VV where V denotes the volume of the 


manifold. In the problem (1), (2) we have A»> 0. 





Received July 4, 1948. 
1We use the usual notations of tensor-calculus (g is the determinant of the covariant metric 
tensor giz). 
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We always assume V together with its boundary (if this exists) to be suffi- 
ciently regular so that those theorems from the theory of eigenvalue-problems 
which are required are valid. 

The aim of this note is to study the analytic continuations in the s-plane of the 
Dirichlet’s series (summation from m = 0 to + or from m = 1 to + @accord- 
ing as Xo > or = 0) 





om(P)dm(Q) 

4 Onh< On\X) 

(4) ) — 

and 

(5) on®(P) 
An’ 


where if V has a boundary, P and Q shall be interior points of V. In the case 
where V is a bounded two-dimensional Euclidean domain, the series (5) was 
first studied by Carleman [1]. Later in the case where V is a bounded Eucli- 
dean domain of arbitrary dimension N it was shown by Minakshisundaram [7] 
by a method different from Carleman’s that (4) is an entire function of s with 
zeros at negative integers and that (5) is a meromorphic function with simple 
pole at s = 4N and zeros at negative integers. The method here developed 
is a generalization of Carleman’s. 

Even though our results are valid to a certain extent under less restrictive 
regularity conditions it is convenient to state them here for an analytic 
manifold V. If A»> 0 (the formulation of the results is only slightly different 
in the case when A» = 0) we find that both the series (4), (5) can be continued arbi- 
trarily far to the left of their abscissas of convergence. The continuation of (4) is 
an entire function with zeros at non-positive integers while (5) represents a func- 
tion analytic except for simple poles at 

s=4N —»,» =0,1,2,.... if N is odd, 
and at 

s = 4N,4N —1,4N —2,...2,1 if N is even. 
The residue at the poles can be determined in terms of the gu. If N is odd the 
function defined by (5) has zeros at non-positive integers and if N is even its values 
in these points can be explicitly determined from the metric tensor of V. By 
Ikehara’s theorem (see [14], p. 44) we obtain as a corollary the relation 


w/2 
yy =P) ~ ae 
a avayn( Y + 1) 





where the sign ~ indicates that the quotient of both the sides tends to 1 when T 
tends to + @ ; see [1]. 


In the case of a closed manifold the series 


(6) 2 Am * 


m=1 


is easily seen to have properties similar to those stated for (5) and by the help 
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of Ikehara’s theorem we obtain immediately the asymptotic distribution of the 
eigenvalues. In the case when V has a boundary our method does not give 
such complete results concerning (6). It is possible by generalizing Carleman’s 
method to deduce the asymptotic eigenvalue-distribution also in this case, but 
it seems as if this could be done more easily by already available methods (see 
[3] and [13]). The analytic continuation of (4), (5) and (6) in the case of a 
sphere was previously studied by Minakshisundaram [8]. 


1. Construction of a parametrix. We introduce normal coordinates in the 
neighbourhood of an inner point P of V. If rpg denotes the geodesic distance 


from P to Q and = differentiation along a geodesic from P, the normal co- 
r 


ordinates (y*) of Q are defined by 
y= ’PQ (= ° 
dr /Q=P 
If @ is a function of r = rpg only and U is an arbitrary function we observe 
that 





> N-—1 d®_ dilogvyg dt dU d@ A 
Frag = o(S + r dr + dr 4) . dr dr +o 
on using the well-known formulae 

ro = git(P)y'y*, 

giz(Q)y* = giz(P)y*, 
where the fundamental tensor is determined with respect to the normal co- 
ordinates and (y*) are the coordinates of the point Q. 

We define in a neighbourhood of P 


a ok. 
2 


a 





H,(P,Q;t) = (Uo + Uit +... Unt”) 


1 
(2./™)" 
by choosing U,(P, Q), » = 0,1, 2,... , independent of ¢ and solutions of the 
differential equations (U_, = 0) 

, aU, , r dlog Ve 7 

dr 2 dr 
The functions U,(P, Q) are uniquely determined by the conditions that they 
shall be finite for P = Q and by the normalizing condition U)(P, P)= 1. It 
is apparent that the choice of the integer m is limited by the regularity of V. 
However, we shall assume V so regular that n can be chosen>%4N —2 (see Sec. 2). 


ot vU, = AU,-1. 





We find 
P)\t 
var.0 ($2) 
o(P, Q) (0) 
and for »y > 0 . ; 
— UoAP,Q) f ren Sn U,-1(P, 11) 
UP, Q) _ req | UP, ll dr py ‘ 


By the help of (7) we find, on account of the definition of U,, that 











wl 
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8 o- 2a. ete” 
For Q = P,t = 0 the singularity of H,(P, Q; t) coincides with the singularity 
of a fundamental solution of the heat-equation 
Au — i 0. 
at 
H,(P, Q; t) is a parametrix of this equation. 
By use of a Laplace-transformation we obtain from H,(P, Q; ¢) the function 


1 wef tHe Be 
K,(P, Q; — &) = Qs y u, |e ad a a. 


The singularity of this function for Q = P coincides with the singularity of a 
fundamental solution of the equation 

(9) Au — tu = 0. 

From (8) it follows that 


AU, r Ho —Ttn 
(10) (A — §)K,(P, Q; — &) = Qvay i ° df dt. 
0 


The function K,(P, Q; —&) is a parametrix of (9). 

The construction of H, (and K,) is analogous to Hadamard’s construction 
of a fundamental solution of the General Wave Equation (see [6]). But by 
comparing our construction with Hadamard’s we see that Hadamard’s proof 
of the convergence of the infinite series he considers cannot be used for the 
infinite series we obtain from H,(P, Q; t) by letting m tend to infinity. The 
same is seen to be true a fortiori for the infinite series obtained from 

K,(P, Q; — &). 
As a parametrix in the large we consider 
rn(P,Q; —§) = nr(rpo)K.(P, Q; —€) 


where npr(r) is a continuous function of r satisfying 


( 1 when r < R 

or(r)= 4 2’ 

| 0 whenr>R, 

and having continuous derivatives of order one and two. In the interval 
4R <r < R the function np(r) can be chosen as a polynomial satisfying in- 
equalities of the form 

dnp < const. dnp const. 

dr|~ R dr? Re - 

R shall be chosen so small that the geodesic sphere round P with radius R is 
contained in the neighbourhood of P where the construction of K,(P, Q; —£) 
is valid. In the case of a closed manifold V we can choose R independent of P. 


"R < 1, < 
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In the case of a manifold with a boundary the choice of R must depend on the 
distance from P to this boundary. 


2. The Green’s function. With the aid of the parametrix I,(P, Q; —&) 
we may express the Green’s function of (9) in the form 
(11) G(P,Q; —§&) = Ta(P, Q; —&) — ya(P, Q; —&). 
The integer m shall be chosen so large, n > 4N — 2, that all singularities of 
G(P, Q; —&) are contained in T,(P, Q; —£) and y,(P, P; —&) is finite. Asa 
function of the point Q the “regular part’’ of the Green’s function, viz., 
vn(P, Q; — &) satisfies the equation 

(4 bed Eoral(P, Q; — £) = (4 — Eo .(P, Q; — €). 
If V has a boundary, G(P, Q; —&) shall satisfy the prescribed boundary con- 
dition (2) or (3). On account of the vanishing of I, and on S the func- 
n 


tion y.(P, Q; —é) satisfies the same boundary condition as the Green's func- 
tion itself. 


Now 7,(P, Q; —&) can be obtained as a solution of the variational problem: 
to search for the minimum of (& is supposed to be real and positive) 
. Ou Ou 
E(u;T,) = sc 2 2Fu) av 
ts: Fa) | (c dx* ax*® lactis: 
where F(Q)=(4 — £)oI'.(P, Q; —£). The admissible functions u are assumed 


to be continuous with piece-wise, continuous, first-order derivatives in the 
open kernel of V. The integral 


[ee S+e)ev 


shall be finite. If the boundary condition u = 0 on S is considered, the admis- 
sible functions shall also satisfy this condition (see [4], p. 482). 
The minimizing function y7,(P,Q; —£) satisfies the equation (see (10)) 


(12) vA(P, P; —§) = Elvan; Tn) -| r,(4l,— §Ta)dV 
V 
from which we observe that on the one hand 
v.(P, Fi —£) < — | r.ar.— tT ,) dV, 
V 
and on the other 
vnl(P, P: — §) = a : | FdV — | r,(4l,.— éT ,) dV 
g Y 
— ; | (aP.— éT.)'dy — | r.(4Pa— él) dV. 
Y 
The first inequality follows from the fact that u = 0 is an admissible function 





~~ —_—— + — 
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in our minimum-problem, the second is obtained by forming a complete square 
under the integral-sign in E(u;I,). Thus an estimate for y,(P, P; —£) for 
large positive values of — can be deduced from the estimate for 





(13) | r,(4T,— éf,) dV 
and 
| (14) ; | (aP.— &',)%4V. 
| By help of the estimations for (13), (14) obtained in the next paragraph we find 
(15) | va(P, P; —&) | < const. — 


where the constant depends on R. 


' 3. Auxiliary estimations. On account of the inequalities 
1 r? r/Jt 
c1t(u+t)2 04 
| |nrl 5 § rr 


2 
, and 
at « 
(16) cC* > |U,| ? < const., (& > &> 0), 
vy=0 
it follows when N > 2 that 
| 
rt 2 ae 
? GF Cc 
(17) |Ta(P,Q;—& | < const. C C™t *dt = const. i 
0 
When N = 2 we use instead of (16) the inequality 


—_# « 
c* ¥ |U,| ? < const. C~™, (E> fo >0;0<a< 4), 
r=0 


| and by the help of well-known properties of the Bessel-function (see [12], pp. 


183, 80, 202) 
é a? nf 
K,(z) = = | e "dr 
2 
0 
| for z tending to 0 and to + © we obtain in this case 
mali 
; (18) IT.(P, Q; —£)| < const.e * log (r+/€). 
For (10) we see that for r < 4R 
AMY Mt Be Fis. <4 
(19) |AT.— é.| < const.e ” | e *t * dt = const. ¢* e * 
0 


| 
| 
provided — 4N+n>-—1. When —4N +2 < —1 we use the method 
; which gave us (17) and (18) and deduce in this way: 
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rvE 

(20) |AT,.— :T | < const. r~¥ +2"? e? when — 4 N+2n< —-1, 

and -- 

(21) =‘ |aTa— &a| < const. e ” log (r/E) when — }N+n = —1. 

In the interval $R < r < R we find inoqualitics for |AT,— éI',| in which the 
RE 


majorizing expressions contain the factor e *. This fact and the inequality 
r > 4R make it possible to give to these inequalities the same forms as (19), 
(20), (21) the constants being now dependent on R. 

By introducing (17) and (19) in the expressions (13) and (14) we find 
N 


| rscar.— tT,)dV| < const. ¢? ‘ 


and 


N 
=—2n-3 
; i” 





(22) oe tr,)*dV | < const. ¢ 





on account of the fact that dV can approximately be substituted by the 
Euclidean volume-element r¥~'dr dQ. Observing (see the beginning of Sec. 2) 
that — N + 2n + 4 isa positive integer it is easily seen that all combinations 
of the inequalities (17), (18) with (19), (20), (21) give the same result (22). 


4. A fundamental formula. Starting from the relation 
G(P, Q; —§) — G(P, Q; —&) = —(& — £0) Jew, ll ; —§) G(0, G; — fo) dVa 
we obtain by repeated application 
? 


(23) G(P,Q; —&) — X (— VE — &0)"G (P, Q; — fo) 


»=0 
= (— 1)°**(& — &)?*" | G(P, 1; —&) G (I, Q; —f) dVn, 
where G (P,Q; —t) = G(P,Q; —&), 
G0 (P,Q; — fy) = | cw M; — fo) G (II, Q; —&) dVu. 


We assume without cetailed discussion that the integral 


(G(P, Q; —&))? dVg 


is finite when? g > [=| and P is an interior point of V. In the case when V 


is a bounded Euclidean domain this follows from the inequality 
£. 
| G(P, Q; —&)| < Ss. 
TPQ 








2 [a] = integer, a — 1< [a]< a. 








a 











Ne OS 
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For other cases we may refer to the works of Giraud [5] and de Rham [11] (for 


the case when V is closed). It follows that for p > 2 [>| the right side of 


(23) can be developed into a convergent series of eigenfunctions, and we have 
’ 


(24) Ta(P,Q; —§) — yalP,Q; —&) — X (— )"(E — &0)"G™ (P, Q; — &) 


v=0 
= (—] p+1 - p+1 om(P)dm(Q) 
— ee et J (m+ £) m+ &)?*# 


m=0 





Since the right side of this relation is finite for Q = P the singular parts of 
r,.(P, Q; —&) must cancel the singular parts of the sum 


? 
(25) 2 (— 1)"(& — &0)’"G (P, Q; —&) . 


What remains after a transition to the limit, Q — P, is the “finite part” of 
r,.(P, Q; —£) for Q = P minus y,(P, P; —£) minus a polynomial in £ of degree 
p contributed by the finite part for Q = P of the sum (25). In order to cal- 
culate the finite part of T,(P, Q; —£&) for Q = P we have to consider the finite 
contributions for r = 0 from 


Ce 
oe —Nis N_, NA» ™ 
| e “t? dt=2? (A) ‘Ky _(rv®), 
2-"— 
0 


where K;(z) denotes the Bessel K-function of order ¢ (see [12], p. 183). If N is 
odd, 2 = N —2v — 2 is odd and ((12], p. 78) 
T 


2 sin rf 





K;(z) = (J_,(2) - I,(z)) ° 


By this formula and by well-known developments of the Bessel functions 
I_,(z) and J,(z) we find that the finite contribution from I,(P, Q; —£) for 
Q = P is, when N is odd 


1 N » 
26 M, * SS a =D -_- — 1 . 
(26) (E, P) aver LT ( a tet ys 


1 
UP, P). 


When NV is even we make use of similar considerations but now based on the 
formula (¢ is a positive integer; see [12], pp. 79 and 80) 


f-1 
f = 2m ft 
rat =1(2) Fan B28)" + -10(8) 
m=0 : 
- (:)" 
, y 2 


wey TE +m + 1) 





z 1 1 
(10g 2 — = v(m + 1) "ooo + m+ ») 
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where 
senor 4 i ate bay 
:".s.. @ 
¥(l) = — y, 


y = Euler’s constant. 
The finite contribution from T,(P, Q; —£) for Q = P obtained in this case is 


¢- mike 
3) M(t, P)-——y ), Piller 2 log 2 
(27) M,(é, )= Gye - (— 1) oa og 
2 


N 1 - N vant 
-ww-¥(5-»)) + ay Br (-gtetaye™. 
N 


rs 
v= 


2 
So we have from (24) 


(28) M,{E, P) = vn(P, P% —) + A,(é, £0, P) 


on?(P) 
m+ &) (m+ £0)???" 





= (— 1s — ¢£.\0+1 
(— ee - "De 
m=0 


where A»(£, fo, P) is a polynomial of degree p in £ with coefficients depending 
on £) and P and where M,(£, P) is equal to (26) or (27) according as N is odd 
or even. In the case when \»> 0 we obtain by performing the transition to 
the limit o— 0 that 


p 
(29) M,(é, P)- xvn(P, Pp: —&) + y APE 


om’(P) 


alfa 1 = +1 —a | 
0 


When Ag= 0 we transfer the first term of the series on the right in (28) to the 
left, take it together with A,(£, &o, P) and note that since the other terms of 
the equation remain finite for > 0, the value of 

o+190°(P) 


A,(é, £o, P) — (— ir (E—£p) geet! 





remains finite and gives an expression of the form 


p 
1 
») A,(P)i — VE 
7=0 
V being the volume of the manifold. We obtain in the case when Ao= 0 the 
formula 








mn 


as 


a a ak Oe 








— —Snn aa 
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? 1 

(30) M(t, P)— v(P, P; 8) + & AUPE - 7 
om’ (P) 
—_— o+1,.9+1 _—— i + a 
as m+ Baer 
m=1 

5. Analytic continuation of (5). We suppose first \»> 0 and multiply both 
sides of (29) by 


1 an a e log | €| ~is(arg E—) 
2mri(— £)* 2m 

and integrate along the following contour in the complex £-plane with a cut 

along the real positive axis. From +© toa (a real, 0 < a < Xo) along the 

lower part of the cut, from a to a along a circle round the origin and then from 

a to +@ along the upper part of the cut. We obtain in this way when the 

real part Rs of s is sufficiently large 








foo] @ ie =] 
on?(P) sin rs f M,C, P) | vx(P, P; —&) 
: 7 — | SX) ; Wie ae 
avy = r ( —— 
? A (P)a’** a'~*e"s* 1 
Ld _— 0 i0(1—8) 
+) eta) Qn | Flae*, P)¢ a 
v=0 0 
where » 
om*(P) 
— (— 1)0+1 gett a. 2s 1 
P(e, P) = (— 1 ¢*! or 
m=0 


The last integral in (31) is an entire function of s vanishing when s equals a 
non-positive integer and so is the expression 





sin 7s y A,(P)a’~**" 


T s—v-l1 


y=0 


On account of our estimation for y,(P, P; —&), (15), the integral 
eo 


| aro — &) dé 


is a regular function in the half-plane Rs > 4N — n — 2. In the case when 
N is odd the first integral in (31) is equal to the following expression (see(26)) 


a 


eo 


o N 
M,(é, P) 1 ( N ) a* ' 
3) - di= - Lj YP, P 
(32) | Yr(-S+e+ —— U,(P, P) 





8 ~ 32 —\N 
é avn a ve 


and in the case when N is even (see (27)) 
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m 7 N 

M,(é, P) 1 (—1)° 

a J eBay y r ‘ ns )( N ) 
a = 2 


,-=7F 
2 v 








(33) -(log a a ar — 2log 2 — y(1) - (X- ’)) UP, P) 
N 2 
ts v9 


Nis¢ 


1 - N a” 
om Te t) nS CLP, FP). 
+ oak ( y +9 +1) N . P) 
3 ae 


2 
All taken together we obtain from (31), (32), (33) the following 
THEOREM. The Dirichlet’s series 








@o 


om*(P) 
(34) | awe. 
m=0 
can be continued to the left of its abscissa of absolute convergence. The function 
¢(s, P) thus obtained can be written 








ed t 
1 U,(P, P) : ; \ 
i(s, P) = —» + R,(s, P) if N is odd, 
(2V*) p r(X -»)(; -*+,) j 
2 2 
and 
N_, 
2 P Pp 
1s, P) = — y UAP, P) + R,(s, P) if N is even, 


—\N 
or LN -Fe) 


where in both cases R,(s, P) is regular in the half-plane Rs > 4N — n — 2. 
When s is equal to a non-positive integer (> 4N — n — 2) we have 


t(s, P) = 0 in the case when N is odd, 





and 1 
is, P) @ T=) vy, PP) che Nisom 
’ ~~ F N ’ 4 . 
(2/n)* F-s 
In the case when Ao= 0 we have to consider the series 
@ 
om*(P) 
m=1 


instead of (34) and to use the relation (30) instead of (29). We obtain essentially f 
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the same theorem but because of the term — 7 in (30) we arrive at the relations 


t0,P) = - when N is odd, 
and 
Uy (P,P) 
0,P) = —2-—, —- - when Ni . 
(0, P) (ava 7 when N is even 


instead of the corresponding relation in the theorem (we denote here by ¢(s, P) the 
analytic continuation of (35)). 

In so far as the Dirichlet’s series (34) or (35) with positive terms can be 
continued to the left of its abscissa of absolute convergence with a simple pole 
at s = $N and residue (we observe that U»(P, P) = 1) 

1 


Qvz)" (*) 
2 
we could apply Ikehara’s theorem (see [18], p. 44) and obtain the following 
asymptotic distribution of the squares of the eigenfunctions as a 








7 
COROLLARY. ») on'(P) ~ . Vv ' 
hm <T (2/n)"T (a + 1) 


6. Analytic continuation of (4). Let P and Q be two different (inner) 
points of V and let R be so chosen that the geodesic spheres round P and Q 
with radius R have no points in common. From the formula (24) we obtain 
when Ao> 0 


? 
(36) Pa(P, Q; —&) — ya(P,Q; —&) + a BP, Q)& 


om(P) m(Q) 


we (am Dot tgots FP Sai’ Om) | 
(— TB Oat Oe! 
m=0 


By choice of R the value of T',,(P, Q; —&) is zero and we have only to estimate 
vn(P, Q; —£€) asa function of §for— + +. Corresponding to (12) we have 


2fya(P, Q; —§) + ra(Q, P; —£)} 
= E(y.(P, 1; —&)+ va(Q, 0; —&) ; T.(P, 0; —§)+ T.(Q, 1; —€)) 
(37) — Elya(P, 1; —&)— ya(Q, 0; —§) ; Ta(P, 1; —&)— T.(Q, 0; —8)) 
+ 2{D(r,(P, 0; —£),1.(Q, 0; —£))+D(r.(Q, 0; —£),T.(P, 1;—£))}, 


where 


D(u,v) = — | wa — tv) dV, 


and where I] is the point of integration. (37) can be deduced in a similar way 
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to (12); see [14]. The last expression on the right of (37) vanishes on account 
of our choice of R. When 

u(Il) = ya(P, 1; —£&) + va(Q, W; —&) 
and 

u(II) = ya(P, 1; —§&) — yva(Q, 0; —6&), 
the expressions 

E(u; T,(P, 0; —£) + 7T.(Q, 0; —&)) 
and 

E(u; rT, (P, i; —&) - r,.(Q, I; — §)) 
attain their minimum values. These minimum values can be estimated as in 
Secs. 2—3 and we find (it is easily seen that we need only the estimation for 
the expression (14)) 

N_sn-3 


(38) | vn(P, Q; —8)+ rn(Q, P; —£)| < const. ¢? 
Interchanging P and Q in (36) and adding we obtain, on using (38) 
N 


TI om(P)m(Q) N_oy— 


? 
3 
(9) (—ayptigt? Yr = 2 AP. e+ OF). 
0 


m= v=0 
From this relation it follows as in Sec. 5 that the Dirichlet’s series 
<1 dn(P)}m(Q) 
(40) > a 
m=0 


which has a finite abscissa of convergence can be continued analytically to the left 
of this line and represents a regular function {(s, P,Q) forRs > 4N — 2n — 3. 
f(s, P, Q) = 0 when s is equal to a non-positive integer (> 4N — 2n — 3). 

When Ao= 0 we consider (40) but with summation from m= 1. Using 
instead of (36) the formula (see Sec. 4, formula (30)) 


? 

1 

P.(P,Q; -8) — (P,Q; -) + BP, OE - 7 
vy=0 


dm(P)om(Q) 


= (—1)*t1+1 
(— tet Do 
m=1 


We find that the analytic continuation of 


(41) y on(P)om(Q) 
m=1 al 


has the same properties as the analytic continuation of (40) with the only exception 
that the function (s, P, Q) represented by (41) is not zero for s = 0 but has the 


1 
value — —. 
V 





of 


ti 


ai 


if 
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Cc 
7. The series > },,’ in the case of a closed manifold. We add a discussion 
m=1 


of the series 
@o 


Lm 


m=1 
in the case of a closed manifold V. In this case R can be chosen independent 


of P and we can make use of the inequality (15) with a constant independent 
of R. Integrating (30) over V we have on account of this inequality 


c 1 
_. et Pies on | ; 
(42) (— tte Yaga = | aeace, Pras 
m=1 V 
t 
ver 
+ X ’ | amy -, +08). 


Using the same method as in Sec. 5 we obtain the 
THEOREM. [If V is closed the Dirichlet’s series 


@ 
Lm 


m=1 


can be continued analytically to the left of its abscissa of convergence and the func- 
tion thus obtained can be written in the form 


. a { U,(P, P) dV 
(2./")* y t )( W ) + R,(s) when N is odd, 
vo T(— — s—-—+y» 








and 
N 


——-1 
a f UMP, P) dV 


(Vay y @ ( y ) + R,(s) when N is even. 


=o T s-—+y 
2 








—-y 


In both cases R,(s) is regular in the half-plane Rs > 4N —n —2. Fors =0 
the value of the analytic continuation is — 1 if N is odd, and 
1 


—. Un(P, P) dV —1 
aoa | ge *) 


Vv 


if N is even. For s equal to a negative integer the analytic continuation is zero 
when N is odd, and when N is even its value is 


ra-s) [, ' 
— Uy (P,P) dV. 
(2+/n)* 1 ;- 
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Ikehara’s theorem gives the asymptotic distribution of the eigenvalues 
N 
VT* 


(2./n)* r(2 + i) 





NQm< T) ~ 


where N(\m< T) denotes the number of eigenvalues < T. 
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ON THE MOTION OF THREE VORTICES 
J. L. SYNGE 


1. Introduction. In a perfect incompressible fluid extending to infinity, 
the determination of the motion of N parallel rectilinear vortex filaments in- 
volves the solution of N non-linear differential equations, each of the first 
order. The method of Kirchhoff' provides certain constants of the motion. 
If we describe the positions of the vortices by their point-traces on a plane 
perpendicular to them, the following facts follow from the theory of Kirchhoff: 


(1.1) The mean centre of the system is fixed. 
(1.2) >’ Kmkn log fmn= const. 
m,n 
(1.3) X Km? = const. 
Here the summations cover the range 1, 2,... N; the prime indicates that 


m =n is omitted; x» are the strengths of the vortices; r,,, is the distance 
between the vortices of strengths x, and x»; rm is the distance of the vortex 
of strength «x, from a fixed point. 

In this paper we shall be concerned solely with the configurations of the 
vortex system, understanding by configuration the geometrical figure formed 
by the vortices, without regard to rigid body displacements of that figure. 
Thus, if a system of three vortices forms a triangle with sides of fixed lengths 
throughout the motion, we say that the configuration is fixed. 

The following theorems, applicable to a system consisting of any number 
of vortices, are obvious from the usual equations of vortex motion, and are 
quoted here for reference. 

THEOREM 1: If, given a configuration, the strengths of all the vortices are 
suddenly reversed, the system retraces the sequence of configurations through which 
it has come. 

THEOREM 2: Given at t = ty a configuration in which all the vortices are col- 
linear, then the configurations at times t = to + r are reflections of one another 
for all values of r. 

THEOREM 3: A system cannot pass through more than two distinct collinear 
configurations ; the times required to pass from one collinear configuration to the 
other are all the same. 

THEOREM 4: Suppose that there are two systems of vortices, S,; and S», each 
consisting of the same number of vortices, and the strengths of the vortices in S, 


Received May 2, 1948. 
1Cf. Sir H. Lamb, Hydrodynamics (Cambridge, 1932), 230; H. Villat, Legons sur la théorie 
des tourbillons (Paris, 1930), 46. 
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being those of S, all multiplied by the same factor K*; suppose further that initially 
the configurations are similar, without reflection, the lengths in S, being those in 
S; all multiplied by the factor L. Then the subsequent configuration of S, after 
time t, is similar, without reflection, to the configuration of S, after time t,, where 
te= t(L*/K*). 

As an immediate consequence of (1.2) and (1.3), we have the following 
result: 


THEOREM 5: [If the strengths of all the vortices have the same sign, their mutual 
distances are bounded above and below for all time, positive and negative. 

No further general results appear to be available, so we turn to special cases. 
The general case can be specialized in a number of ways. We might specialize 
the strengths of vortices, perhaps choosing them all of the same strength, or 
plus and minus one fixed value. On the other hand we might specialize by 
restricting the number of vortices in the system, and this is in fact the speciali- 
zation we shall adopt. 

Since the case of two vortices is trivial, we turn to the case of three vortices, 
without imposing any particular a priori condition on their strengths. This 
is precisely the problem discussed by W. Grébli? over seventy years ago. 
However, he was interested in obtaining formal analytic solutions for the 
motion, and found it necessary at an early stage to specialize the strengths 
of the vortices. He seems to have missed the interesting fact that the motions 
may be classified according to the positive or negative character of the sum of 
the products of the strengths in pairs, xox3+ «xsx1+ «x2. It seems appropriate 
therefore to take up this problem again, concentrating on a qualitative classi- 
fication of all possible motions rather than on the development of analytic 
solutions. The basic equations (2.5) are the same as those of Grdébli, but are 
obtained here in a simpler way. The representation of the motions by trilinear 
coordinates is believed to be new. 


2. The equations of motion and their integrals. Let «1, x2, x; be the strengths 
of the three vortices (i.e. the circulations around them), and R;, R2, R; the 
lengths of the sides of the triangle formed by them, R, being opposite «;, and 
so on, so that, in the notation of (1.2), Ri= 123, etc. In accordance with the 
usual convention, we regard a strength as positive when it gives a counter- 
clockwise circulation. It is convenient to get rid of the factor 2x by defining 
(2.1) Ry= x;/2n, ko= xo/2a, ka= x3/2e. 

It is assumed that none of the three strengths vanishes. 

Consider the rate of increase R’, = dR’,/dt of the side R;. The motions due 
to the vortices k, and k; at its extremities contribute nothing to R’;. One end 
of Ri, viz. ke, has due to k; a velocity of magnitude k,/R; perpendicular to R3;, 
and the other end, viz. ks, has due to k; a velocity of magnitude k,/R, per- 


*Vierteljahrschrift der naturforschenden Gesellschaft in Ziirich, vol. 22 (1877), 37-81, 129-167. 
Grdbli also investigated certain cases of symmetry for N vortices. 
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pendicular to Ry. Let 6;, 62, 6; be the angles of the triangle formed by the 
vortices. Then, on reference to Figure 1, it is seen that 


(2.2) R',; = ek,(R> sin 6;— RR; sin 62), 





Fic. 1 


Rate of growth of a side of the triangle. 


where ¢« = + 1 or —1 according as the circuit of the triangle in the order kykoks 
is positive or negative respectively (counter-clockwise or clockwise). Let A 
denote the area of the triangle, prefixed by a plus or minus sign according as 
the above circuit is positive or negative. Then «A is positive, and 


(2.3) «cA = RR; sin = 4R3R, sin .= 4RiR, sin 63. 
We have also the formula 
(2.4) A= [s(s — R,)(s — R2)(s — R;)]', 


Ss 


4(Rit Rot Rs). 
If we substitute from (2.3) in (2.2) and the two similar equations, we get 
ky RiR’; = 2A(R2?— Rs”), 
(2.5) ko"R2R’, = 2A(Rsv?— Ry), 
ks 'R3R’s = 2A(Ry?*— Rr). 
Adding and integrating, we get 
(2.6) ky R? + ko 'R?? + ks R? = a, 
where a isa constant. If we multiply (2.5) in order by R,*, R,-*, Ry, add, 
and integrate, we get 


(2.7) ky log Ri+ ko log R:+ ks? log R; = b, 
where } is a constant. This is the same as Kirchhoff’s equation (1.2), and 
(2.6) is equivalent to (1.3), but more convenient for our purpose because 
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expressed in terms of the sides of the triangle. The above equations were given 
by Grébli (loc. cit.). 

The differential equations (2.5), with their integrals (2.6) and (2.7), form 
the basis of our work. To these we shall add another equation, obtained by 
differentiating (2.4) and then substituting for R’,, R's, R’s;. In this way we get 


(2.8) A'= F(R, Rs, R;), 
where 


(2.9) f(Ri, Re, Rs) 
=$(2k:Ri'(Ra* — Rs™)][(s — Ri) (s—R2)(s—Rs) + sB(s—R:2)(s—Rs)] 
— s=k,R,"(Rz*— R;*)(s — R:z)(s — Rs). 


Here and later, = indicates summation over a cyclic permutation of suffixes. 


3. Fixed configurations. Let us now seek necessary and sufficient condi- 
tions that the configuration of the three vortices remains fixed, so that the 
motion is a rigid body motion. If the configuration is fixed, then R’; = R’,=R’; 
=0 and so by (2.5) we must have either R, = R, = R; (equilateral configuration), 
or A = 0 (collinear configuration). These are necessary conditions. Any 
equilateral configuration does remain fixed, as was pointed out by Grdbli (loc. 
cit.), and this is a sufficient condition. But A = 0 is not a sufficient condition 
for fixity. At first sight this appears to be in conflict with (2.5). Suppose we 
take for R:, Ro, Rs; any three constant values satisfying one of the equations 
(3.1) Ri= Ro+ Rs, Re= Rs+ Ri, R= Rit Rz, 
such values make A = 0 by (2.4), and hence these values constitute a formal 
solution of (2.5). However, it is a singular solution, and does not in general 
satisfy the full set of equations of vortex motion. In order that the collinear 
configuration may remain fixed, it is further necessary that A’ = 0, or 
(3.2) f(Ri, Re, Rs) = 0, 
where f is as in (2.9). We may sum up as follows: 

THEOREM 6: Necessary and sufficient conditions for a fixed configuration are 
either that the initial configuration be equilateral, or that it be collinear, satisfying 
(3.2). 


4. Variable configurations and the trilinear representation. The values of 
Ri, Re, R; determine a configuration to within a reflection. Thus we might 
discuss changes in configuration by following a representative point in a space 
in which R,, Re, R; are taken as rectangular Cartesian coordinates. Since 
these quantities are essentially positive, we would be concerned only with the 
positive octant. Collision of the representative point with one of the walls 
of this octant would correspond to a collision of two of the vortices. The motion 
of the system would correspond to a curve of intersection of surfaces (2.6) and 
(2.7), the sense in which the curve is described being determined by reference 
to (2.5), with use of the fact that ¢ increases. But the representative point is 


eS eee = 





—_. 
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further restricted since R;, R:, Rs; must always satisfy the triangle inequalities 
(4.1) Ri< Ret Rs, Re& Rs+ Ri, RsX Rit Ro. 

In fact, the planes (3.1) form boundaries in the representative space which the 
representative point is forbidden to cross. If the representative point meets 
one of the planes (3.1), the configuration becomes collinear. Then, by Theorem 
2, the system passes back through the same sequence of configurations but 
with the orientation reversed; the representative point moves back along the 
curve by which it came to the collinear configuration. 


R 














Fic. 2 


The trilinear representation. 


However, there is another and better representation by trilinear coordinates 
in a plane, as shown in Figure 2, and that is the representation which will be 
used in this paper. P,P:2P; is an equilateral triangle of unit height, and x,, 
%2, X3 are trilinear coordinates, i.e. the distances of a general point from the 
sides of the triangle P,P2P;; these values satisfy 
(4.2) Xi+ Xe+ x3= 1. 

Now put 

xi1= Ri(Rit Rot Rs)", 
(4.3) x2= Ro(Ri+ Rt R;)™, 

x3= R3(Rit+ Rot R;)™, 
and so connect the points of the representative plane with the configurations 
of the vortex system. To each configuration of the system there corresponds 
a unique x-point, with one exception: a triple collision (R;= R,= R;= 0) is 
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not represented. On the other hand, to a given x-point there corresponds a 
single infinity of configurations, all similar to one another, together with the 
reflections of those configurations. The centroid E of the triangle (x1= x.= 
xs= 1/3) corresponds to all equilateral configurations. 

Let Q:0203 be the middle points of the sides of the triangle Pi:P2P;. On Q.Q; 
we have x,;= $ and hence x;= x2+ x; or Ri:= R2+ Rs. In fact, Q2Q; corres- 
ponds to the first of (3.1), and the three sides of the triangle 01020; correspond 
to the three planes (3.1) which the representative point is forbidden to cross. 
Since E is certainly permitted, the permitted region is the interior of the 
triangle Q,0:0;. The points Q,0.0; correspond to collisions of the vortices, 
k, and k; colliding at Qi, etc. 

All points on the sides of the triangle Q,0.Q; correspond to collinear con- 
figurations. Since the configuration can change its orientation only by passing 
through a collinear configuration, we may use the two sides of the representative 
plane, all configurations with positive orientation being represented on the 
front of the plane and all configurations with negative orientation on the back. 
The sides 0,020; are then cuts by which the representative point passes from 
one side of the plane to the other. We might in fact throw away all the 
diagram except the triangle 010203, and allow the representative point to pass 
round the edges of this triangle. 

As the system moves, the representative point describes a curve C. To find 


the differential equations of C, we differentiate (4.3) and substitute from (2.5). 
This gives 


(4.4) x’ = KH, x's= KH, x’3= KH, 
where 

(4.5) K = 2AR,?*Rz?Rs"*(ER;)’, 
and 


Hy= — Ryxy(x2?— xs") + xBkixi(x2?— x3"), 
(4.6) H,= — Roxo(x— xy") + XeDk x(x" — x3"), 
Ha= — koxs(xy?— x2") + xsEhkiri(xe?— x;*). 


We check that Hi+ H:+ H;= 0, as of course it must be, by (4.2). 
By (4.4.) we have 


(4.7) —o—oe — oe EA. 


The first two of these equations define a congruence of x-curves, and this con- 
gruence defines the behaviour of the configuration, except for orientation, rate 
of change, and scale. However, orientation is determined by the side of the 
representative plane on which the point lies, and rate of change is given by 
(4.4). As regards scale, if the shape of the configuration is given, its size may 
in general be determined by (2.6) or (2.7), the values of the constants a and b 
being given by the initial configuration. There is, however, one exceptional 
case, and this we shall now discuss. 











ON THE MOTION OF THREE VORTICES 263 


The integrals (2.6) and (2.7) may be written 
(4.8) ki xy + ko x? + ke xe = a(Rit Ret R;)™, 
(4.9) ki log x1+ ka log x2+ ks log xs 

=b —(ky + ko'+ ks) log (Rit Rot R3). 

If a = 0 and 
(4.10) koks+ kakit kika= 0, 
then (Ri+ R:+ R;) disappears from (4.8) and (4.9). In this exceptional case, 
the values of x, x2, x3, a, b fail to determine the values of Ri, R2, Rs. We may 
state the following results. 

THEOREM 7: If the strengths of the vortices do not satisfy (4.10), and b is 
known from an initial configuration, then to each x-point there corresponds by 
(4.9) a unique configuration, except for orientation. 









always singr. x= 
singr. hyperbola pt.at E 


singr. pt.if 
Q, k,+ k, = 0 
Fic. 3 


Singular points. 
(Hyperbola drawn for 2k; = — kz = — k.) 


THEOREM 8: [If the strengths of the vortices satisfy (4.10), then to each x-point 
on the conic 
(4.11) Rokyx+ k3kix2?+ kikox?= 0 
there corresponds a single infinity of similar configurations of both orientations; 
to each x-point lying off the conic (4.11) there corresponds by (4.8) a unique con- 
figuration, except for orientation. 

It is easily seen that, under the condition (4.10), the conic (4.11) is a hyper- 
bola. It passes through the centroid EZ, and meets two sides of the triangle 
0:02:03, each in one point. At E the tangent to (4.11) has the direction 
given by 
(4.12) dx: dx: dxz= ki(ko— kz): Ro(Ra— hi): Ra(Ri— ke). 
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The hyperbola is shown in Figure 3 for the case 
(4.13) 2ki= — ko= — ky. 

It is important to know that no curve C can cut a median of the triangle 00:0; 
in an infinite number of points. To show this, we consider the median P,Q,, 
on which we have 

X= X3= ¥(1 — x). 
By (4.8) and (4.9) we have at an intersection of a curve C with the median P,Q; 
(4.14) ky *xP+ (ko + ks) 3 (1 — x1)?= a/4s*, 
ky log x1+ (ko1+ ks) log (1 — x1) = 6 — (Ri *+ Ra + Rs) log 2s, 
where, as earlier, 2s = R:+ R.+ R;. If we eliminate s, we get an equation 
in x, a, b; for a given curve C the constants a and b are assigned, and this 
equation determines the values of x; corresponding to the intersections of C 


with P,Q;. It is clear that in the range 0 < x:< 1 there can be at most a 
finite number of solutions, and so the result is proved. 


5. Singular points. The most powerful way of studying the congruence 
(4.7) is through its singular points, at which 
(5.1) H,= H.= H;= 0. 

On account of the triangle inequalities (4.1), we are interested only in singular 
points lying inside the triangle 0,020; or on its boundary. Let us first examine 
the points Q;, Q2, Qs, to see if any one of them can be singular. 

At Q; we have x;= 0, x2.= x3;= 4; hence, by (4.6), 

(5.2) Hi= 0, H2= — 4 kot py (ho— ks), Hs = $ bat pe (ka— hs). 
These equations are consistent with (5.1) if, and only if, 

(5.3) ko+ ks= 0. 

When this condition is satisfied, Q, is a singular point. The points Q, and Q; 
may of course be discussed in exactly the same way. 

For all points in the triangle Q,0.0; or on its boundary, other than the 
vertices Q:, Q2, Qs, we have x, x2, x; all different from zero. Then, if we sub- 
stitute in (5.1) from (4.6), we can divide across by these factors, and obtain 
(5.4) x2— x?= ky 8, x?- xr= k28, xv~- x2" = ks, 


= Dkixi (x2? — x3"). 
Addition gives 
(5.5) 6k, = 0. 
Suppose first that @ = 0; then (5.4) give x;= x2= x3= 1/3. Thus the point 
E is a singular point, as is indeed obvious. On the other hand, if (4.10) is 
satisfied, then (5.5) is satisfied with @~ 0. If we multiply (5.4) in order by 
x1", x2", x? and add, we get 
(5.6) Dkr x? = 0 
which is the same equation as (4.11). All singular points (other than Q;, Qo, 
Q3, discussed above) must lie on this conic. Moreover it is easy to see that, 
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if (4.10) is satisfied, then every point on the conic (4.11) or (5.6) is a singular 
point. We have already remarked that this conic is a hyperbola. 

Let us sum up our conclusions about singular points as follows. 

THEOREM 9: The singular points of the congruence (4.7), inside or on the 
triangle Q102Qs, are as follows. If 


(5.7) Rokst+ kakit kiko¥ 0, 

and 

(5.8) Rot ka 0, kat kixX 0, kit koX 0, 

then the only singular point is at E (equilateral configuration). If 
(5.9) Roks+ kakit kike= 0, 


then (5.8) are necessarily true; the singular points make up the hyperbola (4.11), 
which passes through E. If 
(5.10) Rot ks= 0, kat kiX 0, kit koX 0, 
then (5.7) is necessarily true; the only singular points are at E and Q;. Similar 
results hold on permuting suffixes in (5.10). If 
(5.11) ki= — ko= — ky, 
the only singular points are at E, Q2, Qs. Similar results hold on permutation of 
suffixes. 

These results are shown in Figure 3. 


6. Behaviour of representative curves near the point E. To explore the 
curves near the point E, we put 


(6.1) X1= Wit 1/3, x2= yot 1/3, x3s= ys+ 1/3, 
so that 
(6.2) Nt Yot ys= 0. 


Then (4.6) gives, to the first order in y1, yo, Ys, 
Hi= — $ kilys— ys) + Hy Zhily2— ys), 
(6.3) H2= — $ ko(ys— v1) + Hy Thilye— ys), 
H3= — $ ks(yi— yo) + oy Dhilye— ys). 
As in (4.7) we have, as differential equations of the congruence, 


d d és 
(6.4) ata w « 


It is convenient to define 


(6.5) 21= Ya Va, 22= Var V1, 23= Vim J: 

so that, by (6.2), 

(6.6) a - $ (22— 23), = — 4 (23— 21), = — 4 (21— 22). 
From (6.4) we obtain 


(6.7) —s—=—, 
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where 
i I= - £ (H2— Hs) = kote— kts , 


(6.8) I:= — $ (H3— Hi) = ks— kit , 
L;= —- g (Ai- Hy.) = ky2i— kote . 


If we put each fraction in (6.7) equal to ds, we have the equations 


dz 
= = kota— kee, 
dz 

(6.9) lie a 
d. 
= = ky21- Ro&s. 

We have, by (6.5), 

(6.10) 2i+ 22+ 23= 0, 

and so the first two of (6.9) give 
dz 

(6.11) = = kit (ket ks)2s, 
dz 
= = —(kit+t ks)ti— Rate. 

$ 


The solutions are of the form exp (As), where the eigenvalues \ satisfy 
ks— d kot+ kz | 


(6.12) = 0, 
—hi— ks —hs—d | 

or 

(6.13) N= — Thoks. 


Three cases arise: 
Case I: Zkeks> 0; eigenvalues pure imaginary; 
Case II: Zkeks< 0; eigenvalues real, one positive and one negative; 
Case III: Zkeks= 0; eigenvalues both zero. 


7. Case I: kok3+- k3kit+ Riks >0. 


In Case I the curves (6.9) are closed curves, surrounding the point E. How- 
ever, (6.9) is only a linear approximation to the curves C, and it does not follow 
immediately that the curves C are closed. But if a curve C is not closed, then, 
since it cannot intersect itself, it must cut a median P,Q, in an infinite number 
of points. This we have shown earlier to be impossible. Hence all curves C 
near E are in fact closed curves (Figure 4). The sense in which such a curve is 
described depends on the initial orientation of the triangle (cf. (4.4), (4.5)). 
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If we expand the orbit (which, roughly speaking, means bringing two of the 
vortices closer together, since Q:, Q2, Q3 correspond to collisions), we shall reach 
an orbit Co which touches the periphery Q:0.20; at a point corresponding to a 
fixed collinear configuration. This configuration will be approached as a 
limit, not attained in finite time. 


Q 








Q, 


Fic. 4 


Representative curves for Case I: D&ks>0. 


It is interesting to consider here the particular case, ki = k2= ks, which of 
course belongs to Case I. Now the figure is symmetric, and C, will touch all 
three sides of 0:02:03. Thus the system, if started on such a curve, will oscillate 
in infinite time between two fixed collinear configurations, these two configur- 
ations being different. For three equal vortices, the only fixed collinear con- 
figurations are those in which the vortices are equally spaced (Figure 5). Such 
a configuration, if slightly disturbed, will pass in a long time near to one of the 
configurations shown in Figure 6. Equation (2.5) tells us the lengths in Figure 
6 are the same as those in Figure 5. If the representative curve of the dis- 
turbed motion does not meet Q:020; (i.e. if it belongs to the class C of Figure 
4), then all three configurations of Figures 5 and 6 will be approached one after 
another. By symmetry, the representative curve cannot belong to class C; 
or class C,. If it is of class C;, then the motion is an oscillation between a 
collinear configuration adjacent to that shown in Figure 5 and a collinear con- 
figuration adjacent to one of those shown in Figure 6. These oscillations 
between configurations which differ only through interchange of vortices of 
equal strength appear rather interesting. 
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In the general case of unequal strengths, contact will be established first 
with one side of Q:02Q3, as for Cy in Figure 4. When we expand the orbit 
further to C,, we get an oscillation, performed in finite time, between two 
collinear configurations which are actually the same configuration. We may 


think of the return journey as performed on the back of the representative 
plane; it has reversed orientation. 


k, ke ks 





Fic. 5 


Fixed collinear configuration. 
(Ry = k= ks) 


Further expansion gives us C2, which cuts one side of Q:02.0; and touches 
another. Here we have an oscillation between two different collinear con- 


figurations, one of which is a fixed configuration and is not attained in finite 
time. 








Ke ky ki 
ks k, ke 
Fic. 6 


Transforms of configuration of Fic. 5. 


The final stage is C;, representing an oscillation in finite time between two 
different collinear configurations. 


This exhausts the possibilities in Case I. In this case the equilateral con- 
figuration is of course stable for small disturbances. 


8. Case IT: kok3+ k3kit+ Rike< 0. 
Here the eigenvalues are +y, where 


(8.1) u =(— Tkoks)* > 0. 
The solutions of (6.11) are 
(8.2) i= Are" + Bye, 
22> A+ Bue”, 
where 
(8.3) Ailu — k3)- Ax(ko+ ks) = 0, 


Bi(— p — ks) — Balkot+ ks) = 0. 
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As s >, the curve recedes asymptotically in the direction 


(8.4) 2:/Z2= Ai/Az= (ko+ ks)/(u — ks), 
and as s—*+— ~, we have a curve coming in asymptotically from the direction 
(8.5) 2:/22= B,/B2= —(ko+ ks)/(u + Rs). 


These directions may be expressed symmetrically. They correspond to values 
of 21, 22, 3 which make 


dz, dz, dz; 
a da = 21: 22: 23, 
and so, by (6.7), they satisfy 
Ati — Rote + kos = 0, 
(8.6) Rizi + Azo — kys = O, 


— ky2; + Rote + dos 
Zi + 22+ 23= 0. 
If we multiply the first three of these equations in order by k;, ke, ks, and add, 
and then solve with the last of (8.6), we get 
(8.7) 21:22:23; = Mke— k3)+ Skoks— Tkhok; 
: A(Rs— ki) + 3kski— Tkoks 
° Aki ko)+ 3kike— Lkhok;. 

We are to put A = + yu to get the two directions. Figure 7 shows such direc- 
tions (D;, Dz, Dz, D4) and the general nature of the curves near E. 


ll 
= 


< E >D, 





Fic. 7 


Representative curves near E. 
Case II: Dkoks <0. 


The curves which start from E in the directions D,, D2, D;, D, must pass 
out across the periphery Q:10.Q; since they cannot cross nor can they cut a 
median of the triangle an infinite number of times. Similarly all represent- 
ative curves must cross the periphery Q:02Q;. The general nature of the 
pattern is shown in Figure 8. 
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The curves labelled D,, D2, Ds, D4 represent motions in which the configur- 
ation oscillates between the equilateral configuration and a collinear configur- 
ation. The time of approach to E, or recession from it, is infinite. The other 


Q 








Q, 
Fic. 8 


Representative curves for Case Il: Dkoks <0. 


curves represent oscillations between two collinear configurations, not neces- 
sarily distinct. The times involved are finite unless the collinear configur- 
ation involved is a fixed configuration. There are no periodic motions which 
do not include collinear configurations. 

The equilateral configuration is unstable in this case for small disturbances. 


9. Case IIT: Rok3+ k3kit+ kike= 0. 

We have already seen in Theorem 8 that in this case there is a hyperbola 
(4.11) composed of singular points (Hi= H:= H;= 0). If the initial con- 
figuration is represented by a point on this hyperbola, then by (4.4) the repre- 
sentative point remains fixed. Thus the configuration remains fixed in shape. 
To see how it changes its size, we refer to (2.5), in which the right-hand sides 
are now constants. It is clear that the squares of the sides increase or decrease 
linearly with time, remaining fixed in length only if the representative point 
is at E. 

If initially the representative point does not lie on the hyperbola (4.11), then 
both shape and size change. This hyperbola forms a barrier which the repre- 
sentative point cannot cross. Hence the motion consists of an oscillation 
between collinear configurations. 
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ON SURFACE WAVES 
ALEXANDER WEINSTEIN 


1. Introduction. The linearized theory of surface waves leads to several 
mixed boundary value problems which have been investigated by various 
methods. As the physical background of the theory has been repeatedly 
discussed, it will suffice to deal here mainly with the mathematical aspect of 
the question. 

Let D be a finite or infinite domain in the (x, y)-plane and let ¢(x, y) denote 
a function in D satisfying one of the following differential equations 


at) eo 
1.1 oni =—— « QO, 
(1.1) ox? ad oy" 
ao eo 
1.2 — — — Fo = 0 
(1.2) a8 + oy ¢ 


where k? denotes a positive constant. Let m denote the external normal to the 
boundary Cof D. The boundary condition on the part of C corresponding to 
the free surface of the fluid is given by the equation 


(1.3) ot = D9, 
on 
The positive constant p is in some cases an unknown parameter. 
The boundary condition on the part of C corresponding to the rigid part of 
the boundary is given by the equation 


dg 

(1.4) _ 
In the classical theory ¢ denotes (up to a factor depending on the time #) the 
velocity potential in the physical (x, y)-plane. However, in Levi-Civita’s 
theory of plane waves, the independent variables are the velocity potential 
and the stream function. The unknown function ¢ denotes in this case the 
angle which the velocity makes with the horizontal direction. The condition 
(1.3) remains unchanged in form, but (1.4) has to be replaced by the condition 
(1.5) @ = 0. 
In most cases in application the domain D extends to infinity and is bounded 
by straight lines. However in Levi-Civita’s theory of periodic waves, D can 
be mapped conformally on a finite domain such as a circle or circular ring 
without changing the form of the boundary conditions. 

It should be emphasized that (1.3) differs essentially from the boundary 
condition discussed by Fourier in his classical theory of heat conduction. In 
Fourier’s case p is essentially negative, a fact which implies that, for a finite 
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domain D, the corresponding boundary value problem admits only the trivial 
solution ¢ = 0. The situation is, however, different in the case of surface 
waves where the corresponding boundary value problem admits one or even 
several non-trivial solutions for certain positive values of the constant p. The 
important question of the uniqueness of the solution has been overlooked by 
standard treatises on hydrodynamics. Besides its intrinsic mathematical 
interest, a survey of all solutions is of great importance for the following 
reasons: first, the solutions of the linearized problem give a first approximation 
to the exact non-linear theory of a surface wave; second, the superposition of 
two standing waves obtained from two different solutions of the linearized 
problem leads to a travelling wave as required by the theory. 

There are at present three methods of approach to the various boundary 
problems encountered in the theory of surface waves: 

(i) The eigenvalue method. 
(ii) The method of reduction. 
(iii) The method of singular integral equations. 

Some of the problem can at present be discussed only by the first or by the 
third method. However this is not the case for problems discussed up to now 
by the second method alone. It is the purpose of the present paper to show 
that a combination of the method of reduction with the eigenvalue method 
leads to more complete results than the application of the reduction method 
alone. 


2. The eigenvalue method. This method has been developed by A. 
Weinstein [3] in connection with a problem in Levi-Civita’s theory. For 
modification of this method see the papers by G. Hoheisel [4], S. Bochner [5], 
J. L. B. Cooper [6] and A. E. Heins [7]. As an illustration we shall use this 
method for the complete solution of Airy’s Problem which corresponds to the 
hydrodynamical problem of plane waves in water of constant depth: To find 
all harmonic functions in the infinite strip S, - ~- <x<+oe,05ys1 
satisfying the boundary conditions 


(2.1) > os py, fory = 1 
dy 

and 

(2.2) eo = @ for y = 0. 
oy 


Airy’s work contains only a particular solution of this problem which is periodic 
in x and which has been reproduced in all textbooks. 

In order to solve this problem let us consider first the eigenvalue problem 
given by the ordinary differential equation 


(2.3) Y"+ rY = 0, Y= ¥Y(y) 
with the boundary conditions 
(2.4) Y’= py, for y = l, 


(2.5) Y’ = 0, for y = 0. 
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A complete set of eigenfunctions and of corresponding eigenvalues is given by 
the formulas 


(2.6) Yo= cosh apy, Ao= — ag’, 

(2.7) Yn= COS any, An™ Ga’, oe. §. Boca 
where ap is the unique (positive) root of the equation 

(2.8) ao tanh ao= p 

and ai, a2,..., @n, denote the (positive) roots of the equation 

(2.9) a, tan an= p, (m = 1,2,...). 


Turning back to our boundary value problem (1.1), (2.1), (2.2) we develop 
(x, y) for a fixed value of x, into the series 


(2.10) o(x,¥) = x Cn(x) Yn(y). 


This development is possible as ¢ satisfies, for any fixed value of x, the same 
boundary conditions as Y,. The Fourier coefficients c, are given by the 


formulas 
1 


(2.11) Cn(x) = Cy I. o(x, vy) Ya(y)dy 


where 


1 —4 
(2.12) C.= (| Y.'dy) : 
0 


The constant C, is the normalization factor. 
From (2.11) and (1.1) it follows by differentiation that 


1 
c"*,.(x)= — c. | byy Yndy. 
0 


Integrating twice by parts we obtain the formula 
1 
"(x)= _ Caldy Y,.- oY’ nl— c. | oY”, dy. 


0 


The square bracket vanishes in view of the boundary conditions (2.1), (2.2), 
(2.4) and (2.5). By (2.3) we have therefore for c,(x) the differential equation 


(2.13) C'' n(x) — AnCa(x) = O 

which has the following solutions: 

(2.14) Co(x) = do COS ax + bo sin ax 

(2.15) C(x) = ane*** + b,e~*”, (n = 1,2,...). 
On the other hand c,(x) is given by the formula (2.11). 
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Let us assume now that ¢(x, y) satisfies the inequality 
(2.16) I. o(x, y)dy < P4'*!, A>0O 
for |x| +. An application of Schwarz’ inequality to the formulas 


1 
(2.17) Go COS aot + bo sin agx = ca | (x, y) Yo(y)dy 
0 


1 
(2.18) a,e*™" + b,e~*** = c. | (x, y) Y.(y)dy, (n = 1,2,...) 
0 


shows immediately that a,= 5, = 0 for all values of m for which a, is greater 
than A, nm = 0,1,2,.... We have therefore the following result. All solu- 
tions of our boundary value problem satisfying the inequality (2.16) are given 
by the formulas 


(2.19) (x, y) = (ao cos aox + bo sin acer) cosh aoy 
h 
+ ¥ (ane*** + bn,e~*”) cos any 
n=l 


where a and bare arbitrary constants. The exponents a, satisfy the inequality 
(2.20) OS € oct eet... Siat.e << Gest... 

In particular the only bounded solution of the problem is given by the formula 
(2.21) (x, y) = (ado cos aox + D sin ax) cos acy. 


This solution is periodic in x and coincides with the particular solution given 
by Airy. The problem considered in this paragraph cannot be solved by the 
reduction method which will be discussed in the following section. 


3. The reduction method. In this method the unknown function ¢ is re- 
placed by a new unknown (x, y) satisfying the same differential equation 
as ¢@ but vanishing on the boundary of D. The mixed boundary problem is 
reduced to a problem with the classical boundary condition @ = 0. 

T. Boggio [1] was the first to determine all harmonic functions ¢ satisfying 
the condition (1.3) on the boundary of a circle of radius one. Let us put 


(3.1) f(izja=o+ Hy, z=x+ iy = re” 


where y¥ denotes the conjugate function to ¢. Since the real part of ed equals 
Z 


0 . : . , 
r ~ and since the external normal to the circle has the direction of the 
r 


radius r, the harmonic function 


(3.2) =r oe on po 
or 


vanishes by (1.3) on the boundary r = 1. Assuming that ¢ is regular for 
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0 Sr S 1, we see that ¢ vanishes identically and that f sausfies therefore the 
ordinary linear differential equation 
df 


(3.3) .. bf = ta 


where a is a real constant. The integration of this equation shows that a 
regular solution ¢ exists only for p = 1, 2,3,..., in which case ¢ is given by 
the formula 


@ = r?(a cos p§ + B sin p@). 


A more complete analysis of the problem could have been made by the eigen- 
value method, which yields also all solutions ¢ with an isolated singularity 
at the originr = 0. (See Sec. 4.) 

Recently some other interesting mixed boundary value problems corres- 
ponding to waves on sloping beaches have been discussed by Miche [8], H. 
Lewy [9], and J. J. Stoker [10], and others. The method of Lewy and Stoker 
introduces a different reduction procedure. In the following we shall discuss 
as an example one of the problems treated by Stoker and show that a 
combination of the reduction method and of the eigenvalue method yields a 
complete solution of the problem. 


4. A mixed boundary value problem in three-dimensional wave motion. 
Let us consider (see Stoker, loc. cit. [10] paragraph 9) the problem of waves 
in an ocean of infinite depth bounded on one side by a vertical cliff when the 
wave crests are not assumed to be parallel to the shore line. The corresponding 
boundary value problem is the following: 

To find all solutions ¢(x, y) of the differential equation 








ie e 
¢° ¢ 


4.1) 
; 0x? oy 


— ko = 0 


in the domain D, x 2 0, y S 0, satisfying the boundary conditions 


(4.2) o¢ = ¢, for y=0, x>0, 
dy 

(4.3) _. ar forx=0, y <0. 
Ox 


Here k? denotes an arbitrary positive constant. According to Stoker we reduce 
the boundary conditions (4.2), (4.3) to the boundary condition ¢ = 0 by the 
introduction of the function 


0a/fa 
(4.4) a (2 - i) = P(x, y) 
dx \dy 
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which obviously satisfies the differential equation (4.1). It may be written in 
polar coordinates as follows: 
eh 1 a 1 & 

> = 


4.5 — — —- — — PS = 0. 
(4.5) or? r Or t r of 
The boundary condition is 


(4.6) @ = 0, forx = 0, y < Oandy = 0,x>0. 


Instead of prescribing specifically the singularities of @ (see Stoker, loc. cit., 
n. 28, p. 39) we shall determine @ by the eigenvalue method. A subsequent 
integration of the differential equation (4.4) will give us then all possible 
solutions of ¢. 
As in Sec. 2, we consider first the eigenvalue differential problem given by 
the equation 
a6 


4.7 — + r98 = 0 

(4.7) 7 

with the boundary conditions 

(4.8) 8 = 0, for @ = 0 and @ = — 3. 


A complete set of eigenfunctions and of the corresponding eigenvalues is given 
by the formulas 


(4.9) 6,(60) = sin 2n8, X = 4n?, ge 
For a fixed value of r we have for the unknown function @ the expansion 
@ 
(4.10) ®@ = > c,(r) sin 2n0 
n=l 
where 
0 
(4.11) c.(r)= C, | @(r, 0) sin 2n6dé, eS ee eee 
—ir 


where C, is the normalizing factor. From this formula we find by differen- 
tiation with respect to r and by use of (4.5) 





2 —ir /22 
2'.”) + 2¢9) — (e+ my elt) o & | (2 24 sna) sin 2nodé. 
r r r? Jo oF 
The right-hand side in this equation is equal to zero as can easily be seen by 
two successive integrations by parts and by the use of the boundary condition 
(4.6). We have therefore the following differential equation for c,(r) 


2 
(4.12) PP “426. (+ “) cam 0. 
r r? 
The general solution of (4.10) is given in terms of Bessel functions by the 


formula 
(4.13) Cn(r) = Aanlan(kr) + Bant?***Hen™ (ikr) 





~_ —-s An AD Oe [6 


“—™ «— 5, = me see eee =f 86 @& se Oo 
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with arbitrary real constants As, and Bz,. The function J, vanishes for 
r = 0 like r but tends to infinity like e’r~* for r =. The functions #** 
H2, behave like r~" for r tending to zero and tend to zero like e~’r~* at 
infinity. By the same procedure as in Sec. 2 we obtain the following results. 
The solutions @ of (4.5) and (4.6) given by (4.10), can be classified according 
to the behaviour of the integral 


(4.14) L. ®(r, 0)d0. 


The coefficients A», in (4.13) are all equal to zero for any solution # for which 
the integral (4.14) is o(e*’r~) at infinity. The coefficients B;, vanish for n > h 
for all solutions @ for which the integral (4.14) is o(r~**) at the origin. The 
only solution ® which is regular everywhere is ® = 0. 

By taking = 0 and ® = iH,” (ikr) sin 26 and by integrating the corres- 
ponding differential equations (4.4) for ¢, Stoker obtains two standing waves 
which can be combined into a travelling wave. One of these standing waves 
has a logarithmic singularity at the origin. From the results of the present 
paper we see the presence of a singularity is an unavoidable consequence of 
the linearized theory of surface waves. The contradiction of the original 
assumption of small amplitudes is somewhat mitigated by taking the 
solution with the weakest singularity at the origin. From the mathematical 
viewpoint, however, there is no reason to introduce any limitations on the 
behaviour of the solutions. 


5. The method of singular integral equations. We conclude with a few 

remarks about this method which has been applied to the case when the 
domain D is a parallel strip, as in Sec. 2. Let us replace in Airy’s problem 
the condition (2.2) by the condition 
(5.1) @ = 0, fory = 0. 
Under certain restrictive assumptions on the behaviour of ¢ at infinity the 
modified problem can be reduced to a Picard integral equation [2]. However, 
as has been mentioned in Sec. 2, the eigenvalue method gives the solution of 
the same problem under less restrictive conditions. The situation is, however, 
different in the dock problem in a channel of finite depth, which is obtained by 
imposing the condition (2.1) for y = 1, x > 0 and the condition (2.2) on the 
remaining part of the boundary of the strip. This problem, which seems at 
present inaccessible by any other method, has been solved by A. E. Heins [11] 
by a reduction of the problem to a Wiener-Hopf equation. The assumptions 
which are required in order that this problem be formulated as a Wiener-Hopf 
integral equation are discussed in paragraph 9 of the paper by Heins. Unique- 
ness is studied in relation to the Wiener-Hopf integral equation to be solved. 
This integral equation is equivalent to the original boundary value problem 
subject to the conditions mentioned above. The general uniqueness theorem 
under less restrictive conditions, has not been discussed yet for the dock 
problem. 
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ANGULAR MEASURE AND INTEGRAL CURVATURE 
HERBERT BUSEMANN 


Tue Gauss-Bonnet Theorem leads through well known arguments to the fact 
that the integral curvature!’ of a two-dimensional closed orientable manifold M 
of genus p equals 4x(1 — p). This implies, for instance, that the Gauss curva- 
ture' K can neither be everywhere positive nor everywhere negative, if M is 
homeomorphic to a torus. 

The relations between the sign of K and the topological structure of M have 
been the subject of many investigations. Those of Cohn-Vossen [4, 5] are par- 
ticularly interesting, because they are not restricted to closed manifolds. 

Hadamard [6] showed that the condition K < 0 determines to a great extent 
the shape of the geodesics (on closed or open manifolds). The already men- 
tioned papers of Cohn-Vossen show also how the condition K > 0 influences 
the behaviour of the geodesics. 

All these investigations rest on the Gauss-Bonnet Theorem, which states in 
its most primitive form that the integral curvature of a geodesic triangle equals 
the spherical excess of the triangle. Thus they depend ultimately on the concept 
of angular measure. This concept is in turn derived from the local, that is the 
Euclidean geometry, where it means amount of rotation. 

The Minkowskian geometry is the local geometry of non-Riemannian 
metric spaces. It does not permit general rotations. If the distance is sym- 
metric, which will always be assumed here, the Minkowskian geometry permits 
reflection in a point, which in the Euclidean case is equivalent to rotation 
through x. Therefore no particular angular measure can be entirely natural 
in Minkowskian geometry. This is evidenced by the innumerable attempts 
to define such a measure, none of which found general acceptance. 

Of course, it is generally agreed that angular measure must be additive for 
angles with the same vertex. In view of our previous observation, it is natural 
to add the requirement that straight angles have measure x. It will be shown 
here that any angular measure with these two properties permits us to establish for 
general spaces most of the above quoted results on Riemann spaces, provided we 
interpret conditions like K > 0 on M to mean that every non-degenerate small 
geodesic triangle on M has positive spherical excess. For some results it is 
necessary to add a condition, which is always satisfied by the ordinary angles 


Received October 8, 1948. 

1The English expression ‘‘total curvature” corresponds to the German ‘“‘Gauss’sche Kriim- 
mung,” whereas the German expression “Totalkriimmung” is “integral curvature”’ in English. 
In order not to confuse the reader interested in the original literature, the present paper avoids 
“total” altogether by using “Gauss curvature” corresponding to the German, and “integral 
curvature” corresponding to the English custom. 
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in Riemann spaces, and which states essentially that, ina uniform way, an angle 
cannot be nearly straight without having a measure close to rf. 

The main point of the present paper is the tenet that angular measure in 
Finsler spaces is—contrary to the prevailing views—a very fruitful concept, and 
that it becomes unnatural and barren only through insistence on particular mea- 
sures. 

The extent of the material in the Riemannian case precluded its full dis- 
cussion here. Except for a glance at the connection of excess with the theory 
of parallels (Sec. 2) and the topological structure of compact manifolds (Sec. 3) 
the paper concentrates on the work of Cohn-Vossen, whose arguments are 
partly reproduced here. 

Hadamard’s results are only briefly touched because the author showed 
recently in [3], although not in connection with angular measure, that they do 
not depend on the Riemannian character of the metric. 


1. Angles in systems of plane curves. In the Euclidean plane E with the 
(Euclidean) distance xy let S be a system of curves with the following pro- 
perties: 

I. Each curve is an open Jordan curve, that is, it has a representation q(t), 
—2<t<© where q(t) is continuous and q(t;) ¥ g(t2) when t, # te. 

Il. g(t)qg(@0) +© when |t| +@. 

III. Amy two distinct points of E lie on exactly one curve of S. 

The curve in S (S-curve) determined by the two points a, } will be denoted 
by g(a, 6). On g(a, 5) the points a, b bound an arc t(a, 6). The symbol (acb) 
means that a, b, c are three different points and that c lies on t(a, 6). We also 
put t(a, a)= a. 

The S-curves satisfy all the axioms of order and connection of Hilbert, in 
particular the axiom of Pasch.? In addition, a,— a and b,— b implies t(a,, b,)— 
t(a, b) and, if a # 5, also g(a,, b,)—+ g(a, 6). The arrow indicates here that 
t(a, b) or g(a, 6) is Hausdorff’s closed limit of the sets t(a,, b,) or g(a,, b,).2 A 
point p of an S-curve g divides g into two (closed) rays t, t2, which we call 
opposite. 

If tr; and tm are two different rays with the same origin , then 1, U te 
divides E into two (closed) domains D, and Dz. The sets of all rays with origin 
p in D; and D, respectively are the two angles with legs t, and t2. They are 
called straight if r, and rz are opposite. Otherwise exactly one of the domains 
is S-convex* and we call the corresponding angle the convex angle (r:, t2)"™, 
and the other the concave angle (r:, t2)“”. It is convenient to complete this 
definition by letting (r,, r:)"** mean the set consisting of r; alone and (1, r;)“*” 
the set of all rays with origin p. Ifa, b, c are three points not on one S-curve, 
then Z abc means the convex angle whose legs are the rays from } through 
a and ¢. 


*Proofs are found in [1, Sec. III.3]. 
%A set X is S-convex if a, b €X implies t(a, })C X. Compare [1, Sec. III.3]. 
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We now assume that an angular measure |D| has been defined for the angles 
D in S with the following properties: 

1) |D| 2 0. 

2) |\D| = x if and only if D is straight. 

3) If D; and Dz are two angles with a common leg but with no other common 
ray, then |\D,\UDs3| = |D,| + |D,I. 

We say that the angle D, tends to the angle D, if the legs of D, tend to the 
legs of D, and if r,eD, and r,— r implies reD. We call the angular metric 
continuous if 

4) D,— D implies \D,| — \D). 

Some consequences of 1), 2), 3) are 
a) \D| = 0 if and only if the legs of D coincide. 

For if the legs r:, tz of D coincide, let D’ denote one of the two straight angles 
with t1= tz as one leg. Then DU D’= D’, therefore by 2) and 3) 

«= |DUD"' = |D| + |D'| = |D| +x, 

so that |D| = 0. Conversely let |D| =. Its legs r, and r, cannot be opposite 
by 2). Denote the opposite ray tor, by r;._ If r; and rz did not coincide and 
D = (t,, t2)” then by 1), 2),3) |D| = + |(ts, r2)""|> 9. If D =(ti, t2)"™ 
then # = |D| + |(r2, ts)"*| = |(r2, t3)"| although (te, ts) 
b) Convex angles have measure less than # and conversely. 

Concave angles have measure greater than x and conversely. 
c) Vertical angles are equal. 
d) The sum of the measures of the angles in a triangle abc (set bounded by 
three segments t(a, 5), t(b, c), t(c, a), where a, b, c are not on one S-curve) is 
positive and less than 3z. 
e) If the angular metric is continuous and the points a,, b,, c, are not on an S-curve 
and tend to a point p, then the sum of the angular measures in the triangle a,b,c, 
tends to x. 

A proof follows immediately from the observation that g(a,, ,), g(b,, c,) and 
a(c,, @,) may be assumed to converge. Then Z b,a,c, and the vertical angles 
to Z a,b,c, and Z a,c,b, tend to three angles whose union is a straight angle. 

Some of the preceding remarks extend in the usual way to degenerate 
triangles and will be used for such triangles. 

The excess ¢(abc) of the triangle abc is defined as 
(1) e(abc) = \abc| + \bca| + |\cab| — x, 
where |abc| = | Z abc| . 

Degenerate triangles have excess 0. If the triangle abc is decomposed (sim- 
plicially by S-curves) into the triangles a,b,c, then 

e(abc) = Le(a,b,c,). 

If a;,..., @, are the vertices of a simple closed polygon P with sides 

t(a;, @;41) and a; is the measure of the angle at a; measured inside the closed‘ 





vex 





is not straight. 


‘In this paper domains bounded by geodesic polygons are always understood to be closed. 
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domain G bounded by P then for any simplicial subdivision of G (by S-curves 
is always understood) into triangles a,b,c, 


(2) e(a,b,c,) = 2 — Z(e — as) = La; — (mn — 2)z. 

Let r; and rz be two opposite rays with origin p determining the two straight 
angles D, and Dz. If a,er;, a; pand geD,—(t, Ure) then a ray r, with origin p 
and through a point xet(a;, g) U t(g, @2) traverses monotonically all rays in D, as 
x traverses t(a:, g)+ t(q, a2). Therefore |(r:, r.)"| = ¢(r,) is a strictly in- 
creasing function with ¢(r;) = 0, $(t2)= 7. The values of (rz) = |(t:, t2)*"| 


for r,eD; are determined by 2). If D is any angle with vertex p which does not 
contain 7; and has legs r’, r’’, then 


(3) D| = |¢(e’) — o(e”)|. 
If D contains r,, then 
(3’) |D| = 2x — |e(r’) — o(r’”’)|. 
Conversely, if in D, any strictly increasing function ¢(r,z) with ¢(r:) = 0, 


¢(t2) = x is given, and ¢(r,) is determined in D, to satisfy 2), then (3) and (3’) 
determine an angular measure at p which satisfies 1), 2), 3). 





2. Excess and parallels. The present section is concerned with the relation 
of the angular metric in a system S to the theory of parallels. It will not be 
needed later on but will elucidate the meaning of an angular metric. 

If g* is an oriented S-curve and x traverses g* in the positive sense, then the 
line g(p, x) converges for any fixed point p toa line a. If g*(a, 6) denotes 
generally the line g(a, 6) with the orientation in which } follows a then g*(p, x) 
tends to an orientation a* of a. We call a(a*) the (oriented) asymptote to g* 
through » (for a proof of this and the next statements see [1, Sec. III.3]). The 
line a does not intersect g. The asymptote to g* through any point gea is again 
a. But in general g* is not an asymptote to a*, for an example see [1, Sec. III.5]. 

Let the parallel axiom hold, that means, through a given point p not on a 
given line g there is exactly one line § which does not intersect g. If we deter- 
mine angular measure at one point as at the end of the preceding section but 
with a continuous ¢, and define measure for an arbitrary angle as equal to the 
corresponding angle at p with legs parallel to the given angle, then condition 
4) is also satisfied and the excess of any triangle is 0. 

(4) A system S in which the parallel axiom holds possesses continuous angular 
metrics with excess 0. 

However, it is not true that zero excess implies the parallel axiom nor does 
the parallel axiom imply that every continuous angular metric has excess 0. 
The only statement which holds without further conditions on the angular 
metric is the following: 

(5) If the excess is non-positive and the angular metric is continuous then the 
parallel axiom implies zero excess. 


Let abc be a non-degenerate triangle. If (abx) and x traverses g* = g*(a, 5) 














——— + 





ANGULAR MEASURE AND INTEGRAL CURVATURE 283 


then |bac| + \acx| < and gt(c, x) tends to the asymptote §* to g*. If y 
follows c on §* then |acx| — \acy| = \acb| + |bcy| so that 


(6) \bac| + \acy| < r. 








aii aa 


If (ycy’) and (bax’) then g*(c, y’) is, because of the parallel axiom, the asymptote 
to g*(b, a) through c. As before we see 


(6’) |x’ac| + |acy’| < + 

and since |x’ac| = x — |bac|, |acy’| = x — |acy| it follows from (6) and (6’) 
that |bac| + |acy| = |x’ac| + \acy’| = so that |bac| = |acy’|. For the same 
reason |abc| = |bcy| and since |acy’| + |acb| + lbcy| = x the theorem is 


proved. 

It is clear that this argument and the additivity of the excess yield the fol- 
lowing more general fact: 

(7) If the excess is non-positive and the metric is continuous and there is only 
one line } through a point p not on q which does not intersect g, then any 
triangle, whose vertices are in the closed strip bounded by g and 6, has excess 0. 

The arguments in the above proof could be reversed if |cxb| +0. The 
following examples will show that further progress is impossible without this 
property. A hemisphere H without the bounding great circle can be mapped 
on the Euclidean plane E in such a way that the arcs of great circles in H go 
into the Euclidean straight lines in E. If we assign toan angle in E the 
measure of the corresponding spherical angle (in the usual sense), then the 
excess in any triangle of E is positive in spite of the parallel axiom. With 
this measure the same holds for the straight line pieces in the interior of a 
circle in E. On the other hand the Euclidean angles in E may be used as 
angular measure for those same pieces. This means that both positive and 
zero excess are compatible with the hyperbolic parallel axiom. 

We call an angular metric in a system S of curves complete if it is continuous 
and | pxq| — 0 whenever x traverses a ray with origin p from p toward @. 
(8) In a complete angular metric the excess cannot always be positive. 

With the same notation as above positive excess would yield ¢(x’cx) > (abc) 
> x. Because of |cxx’| + 0 and |cx’x| — 0 it would follow that for x and x’ 
which are sufficiently far away |x’cx| > « which is impossible. 
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(9) In a complete angular metric zero excess implies the parallel axiom. 

For then |cax| + |acx| + |axc| = x and |axc| +0. Moreover t(c, x) tends 
to a ray t which lies on the asymptote r* through c to g*. If zer, z ¥ c then 
|acx| —» |acz|, hence |cax| + \acs| = x. Similarly t(c, x’) tends to ray r’ on the 
asymptote through c to g*(b, a) and if s’e’, 2’ c, then |cax’| + |acz’| = x. 
It follows from |cax| + |cax’| = x that |acz| + |acz’| = x, so that rand r’ are 
opposite rays, q.e.d. 

(10) In a complete angular metric with non-positive excess asymptotes are sym- 
metric. 

If r* is an asymptote to g* and g* were not an asymptote to r* then the 
asymptote to r* through a point } of g* would be a line f* different from g* 
(see Figure). If cer* and u follows 6 on f* then g(c, u) intersects g* by the 
definition of asymptotes in a point x with (cux). Because the excess in bux is 
non-positive 

leub| > |ubx| + |bxu| > |ubx| > 0 


but |cub| + 0 when u traverses f+ in the positive sense. 

Example 1) in [1, Sec. III. 5] yields, with the ordinary Euclidean angles, a 
complete angular metric and non-symmetric asymptotes which shows that (10) 
would not hold without the assumption that the excess is non-positive. 

(9) and (10) are first examples of statements which connect conditions or 
the excess with topological properties (in this case of the system 5S). 


3. Angular measure for curve systems on two-dimensional manifolds. The 
word surface will be used here to denote a connected two-dimensional topo- 
logical manifold. 

As in Sec. 2 for the plane we consider on a given surface M a system S of 
curves with the topological properties of geodesics. The existence of such a 
system is guaranteed by the following two conditions. 

1) Every point p of M has a neighbourhood U(p) homeomor phic to the plane, in 
which a system S, of curves is distinguished with the properties 1, 11, 111 of Sec. 1. 

2) If a, b, c ie in U(p) (\ U(qg) then (abc) holds with respect to S, if and only if 
tt holds with respect to S,. 

By 2) a segment t (a, b) in S, is also a segment in S,. Therefore the notation 
t(a, b) can be used without reference to a definite system S, as long as a and b 
both lie in some U(). 

The concept of a geodesic will actually not be used in the sequel. But since 
1) and 2) are derived from this concept, we mention that an S-geodesic is to be 
defined as a continuous curve x(t), — © < t <@ with the following property: 
if tp is given and x(to) «U(p) then a suitable subarc t,;< t < ty with th< to < fy 
of x(t) represents a curve in S,. The existence of geodesics can be established 
by the procedure of [2, Sec. II.5]. 

If a eU(p) (1) U(g) then a ray rt, with origin a in S, will, in general, not be 
a ray in S,, but by 2) the ray r, either contains, or is contained in, a ray t, with 
origin a of S,, which is uniquely determined by rp. 














ANGULAR MEASURE AND INTEGRAL CURVATURE 285 


If r,', t,” are two rays with origin a in S, and r,', r7 are the corresponding 
rays in S,, then 2) clearly implies that r,' and r,? are opposite if and only if 
r,' and r,? are. Also, if r,' and r,* are not opposite, then a ray in (r,', r,*)"™ 
or (r,', tp”) corresponds to a ray in (r,', r,*)"™ or (r,', r,*)™” respectively. 

These facts lead to the following formal definition of a ray t with origin a in 
S: tis a set of rays in the local curve systems with these properties: 

a) t contains exactly one ray of every S, for which aeU() and no ray of any 
other S,. 
b) If t,er and r,er then either r, Dr, or tpC fz. 

The meaning of angles in S, of symbols like (r:, r2)"*, and of convergence of 
angles is now obvious. 

An angular measure for the angles in S is then characterized by the properties 
1), 2),3) of Sec. 1 and a continuous angular measure by the additional property 
4). Through the natural one-to-one correspondence between the angles with 
vertex a in S and the angles with vertex a in S,, U(p) Da, an angular measure 
in S induces an angular measure in S,. 

Whenever the word triangle is used it is understood that its vertices lie in 
one U(p). The excess of a triangle is still defined by (1). A geodesic polygon 


P on M isa curve of the form U t(@i, 2:41), 2: G41. Some of the angles of P 


i=l 
may be straight, that is the segments t(a;_,, a;) and t(a;, a;4:) may belong to 
opposite rays with origin a;. 

We then call a; an improper vertex of P, otherwise a proper vertex. If all 
vertices of P are improper then P is a geodesic arc. If in addition a,;= a,4, 
and the angle at a, is straight, we call P a closed geodesic. 

Let G be a compact domain of finite genus on M which is bounded by n 
simple closed mutually non-intersecting geodesic polygons. If G is simplicially 
divided into triangles, then the number of vertices minus the number of sides 
plus the number of triangles is an integer X(G) which depends only on G and 
not on the choice of simplicial division. According to the terminology pre- 
vailing in topology, X(G) is the negative Euler characteristic of G.5 

Any simplicial division of G into triangles a,b,c, satisfies the following 
fundamental relation 


(11) Ce(a,b,c,) = 2nX(G) — X(w — ai), 


where a; are the angular measures of the angles at the vertices of the boundary 
of G measured in G.* It is immaterial whether a; traverses the angles at all or 
only the proper vertices. 

If M is compact and has finite genus p we find that for any simplicial decom- 
position of M into triangles a,b,c, 


5Compare Kerekjarto [8] and Seifert-Threlfall (9). Cohn-Vossen [4] calls X(G) (and not 
—X(G)) the characteristic of G. 

‘A modification of the topological proof for (11) which is adapted to the present conditions 
is found in [4, p. 120]. 
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4x(1 — p) if M is orientable, 
(12) Leb td = 2eX(M) = re — ») if M is not orientable. 
The number 2e(a,b,c,) in (11) or (12) which is independent of the simplicial 
division is called the integral curvature C(G) of Gor C(M) of M. Wesay that 
M or a domain G on M has positive, negative, non-positive, non-negative, or 
zero curvature if for every non-degenerate triangle abc in M or G 


e(abc)> 0, <0, £ 0, 2 0, or = O respectively. 


If the curvature of a twodimensional Riemann space R is non-negative, non- 
positive, zero, positive, or negative in the usual sense then R has the same 
property in the present sense. The converse is true in the first three cases, 
but not always in the last two. If the Gauss curvature of R is positive (nega- 
tive) except on some curves or isolated points, R has still positive (negative) 
curvature in the present sense. The existence part of the following theorem 
follows therefore from well-known facts regarding Riemann spaces, the re- 
mainder is a consequence of (12). 

(13) THEOREM. A compact surface M can be provided with a system S of geo- 
desics and an angular measure such that curvature is: 
a) non-negative, if and only if M is homeomorphic to the sphere, torus, one- 
sided torus (also called Klein-Bottle), or the projective plane. 
b) non-positive, if and only if M is not homeomorphic to the sphere or the 
projective plane. 
positive, if and only if M is homeomorphic to the sphere or the projective 
plane. 
negative, if and only if M is not homeomorphic to the sphere, torus, one- 
sided torus or the projective plane. 
A torus or one-sided torus with non-positive or non-negative curvature 
has curvature 0. 


~— 


c 


d 


~a 


4. Two dimensional metric manifolds. No statements which approach (13) 
in completeness seem to be possible for non-compact surfaces unless the curves 
in S are really geodesics in the metric sense, and not only curves with the 
topological properties of geodesics. That M is a space with metric geodesics 
is expressed by the following conditions: 

A. M is a metric space with distance xy. 

B. M is finitely compact, or a bounded sequence has an accumulation point. 

The fact that the three points a, b, c are different and satisfy the relation 
ab + bc = ac will be written as (abc). 

C. M is convex, that is, for any two distinct points a, c a point b with (abc) 

exists. 

D. Prolongation is locally possible, or for every point p there is a p(p)> 0 
such that for any two different points a;, a2 with a;p < p(p) a point d 
with (a;d@2d) exists. 

E. Prolongation is unique, or, if (@:a2d’), (@,:a2d") and aed’ = a2d” then 
d’= d". 
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These axioms guarantee the existence of geodesics (compare [2]). In the 
present case we add 

F. M has dimension 2 (in the sense of Menger-Urysohn). 

It can be proved that M is a connected topological manifold or a surface 
(for this and the following statements see [1, Sec. 1.4]). A space which satisfies 
Axioms A to F will be called a G-surface. 

A metric segment is an isometric map of a Euclidean segment. If U(p) is 
the interior of a sufficiently small geodesic triangle on M, then the open metric 
segments in U(p) with endpoints on the boundary of U(p) form a curve system 
S, with properties 1) and 2) of Sec. 3. 

Since any two points of M can be connected by a metric segment, only those 
metric segments are segments in the previous sense, which lie entirely in one 
U(p). But since every metric segment can be divided into a finite number of 
metric segments each of which lies in one U(p), and the angles at the points 
of division are straight, the distinction between the two kinds of segments 
turns out to be immaterial and will therefore be dropped. 

If M has finite connectivity it can be represented topologically as a compact 
manifold M of finite genus which has been punctured at a finite number of 
points 2;,..., 2%. Let P be a simple closed geodesic polygon which bounds 
on M a simply connected closed domain T which contains exactly one z;, say 
Zi). Because of B the set T — 2;, appears on M asa set which looks like a half 
cylinder and extends to @. Wecall T a tube (Fluchtgebiet in the terminology 
of Cohn-Vossen [4]). 

The tubes are the new feature of non-compact M as compared to compact 
surfaces. The study of non-compact M must therefore be based on the prop- 
erties of tubes. The remainder of this section investigates tubes. 

With the above notation, consider on T the class C(u), u > 0, of all curves 
C which are homotopic to P on T and have distance at most u from P. Whether 
this distance is measured on M or on T is immaterial. For if measured on M 
then a segment connecting a point of P to a point of C exists whose length 
equals the distance of P and Con M. This segment cannot contain a second 
point of P and lies therefore entirely in T. 

C(u) contains curves of finite length (for instance P). Since 7, considered 
as space, satisfies B and every member of C(u) contains a point whose distance 
from P is at most u, there is a shortest curve R(u) in C(u) (for a proof compare 
[1, p. 10] and [2, p. 234]). The length A(u) of R(u) is obviously a non-increasing 
function of u and the triangle inequality yields easily that \(u) is continuous 
(see [4, §16]). We represent R(u) with the arc length ¢ as parameter in the 
form x(t), 0 < t < A(u), x(O) = x(A(u)). Notice first 
(14) If x(to) is not a vertex of P and has either distance greater than u from P or 
is not the only point of R(u) whose distance from P is at most u then the subarc? 
to— 6 St < tot 6 of x(t) is a segment for sufficiently small 6 > 0. 


?This inequality is to be replaced by the two inequalities 0 < ¢ £6 and A(u)— 5 Ct CA(u) 
if to = O or to= A(u). 
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For otherwise the subarc to— 5 < ¢ < to+ 6 can be replaced by a segment 
with the same endpoints. If 4 > 0 is small enough, the new curve R’ will still 
lie in T, even when x(t) lies on P, but is not a vertex of P. Moreover R’ will 
still be homotopic to P and have distance < u from P. But the length of R’ 
would be less than A(u) which contradicts the definition of R(x). 

(14) implies that R(u) is a geodesic polygon. Moreover, if R(u) contains 
points with distance < u from P then R(u) contains infinitely many such 
points and none of them can be a vertex of R(u). Therefore we see 
(15) R(u) is either a closed geodesic, or all its vertices are vertices of P, or R(u) 
has exactly one vertex and its distance from P equals u, whereas all other points 
of R(u) have greater distance from P than u. 

We show next that R(u) is a Jordan curve. Since R(u) is homotopic to P 
and T is homeomorphic to a halfcylinder, R(w) must contain subpolygon R’ 
which is a Jordan curve and homotopic to P. If R # R’ then R’ cannot have 
distance < u from P, otherwise R’ would belong to C(u) and be shorter than 
R(u). Therefore R(u)— R’ contains a point r with distance < u from P; r 
may be chosen as x(0). Then R’ is a subarc of x(t) of the form 0 <a {< t < 
B < A(u) with x(a) = x(8). ThearcsO < ¢ < aand®B < t < A(u) of x(t) must 
have the same length, otherwise replacing the longer by the shorter would 
yield a curve in C(u) with smaller length than A(w). 

Replacing the arc 8 < ¢ < (u) by the arc 0 < ¢ < a, that is defining y(t) = 
x(t) for 0 [ t < B and y(t + B)= x(a — 2) for Of t K a = A(u)— 8B, yields 
again a curve R* in C(u) of length A(u). Statement (14) would then apply to 
R*, hence for small 6 > 0 the arcs a — 6 <¢t Ka +6 and B —6 Lt< B+ 
would be segments. By construction the arcs a — 6 <¢ Ka and B Lt< pti 
coincide. The uniqueness of the prolongation E would imply that the arcs 
a <t{a+éand B— 5 < t & B also coincide, but then R’ would not be a 
Jordan curve. 

Since R(u) is a simple closed geodesic polygon homotopic to P it bounds a 
subtube 7(u) of T. If avertexr of R(u) is a vertex of P, then the angle of R(u) 
at r measured in 7(u) cannot be convex, otherwise R(u) could, because of 
u > 0, be shortened without violating the conditions for belonging to C(x). 

Finally it will be proved that in case R(u) has exactly one vertex g with 
distance u from P, the angle at g measured in T(u) must be convex. Let t be 
a segment of length u connecting g toa pointdon P. Then t cannot contain 
other points of either P or R(u) because the distance of R(u) from P would 
then be smaller than u, contrary to (15). If the angle D at gin T(u) were con- 
cave, let ¢:, ¢: be points on the legs of D and close to g. Then the interior J 
of the triangle gcicz would lie outside of T(u). Also, J \U T(u) contains a 
neighbourhood of g. The segment t connects d to g without entering T(x). 
It must therefore cross t(c:, c2) at a point g’ and q’ has distance u’<u from P. 
If then the arc t(c:, g) Ut(g, cz) of R(u) is replaced by t(c:, cz), the length 
decreases so that A(u’)< A(u), which is impossible. 

Thus we have proved the Theorem of Cohn-Vossen: 

(16) R(u) is a simple closed geodesic polygon. It is either a closed geodesic, or all 





——— or 
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its proper vertices are also vertices of P and the corresponding angles measured in 
T(u) are concave, or R(u) has exactly one proper vertex q, which is the only point 
on R(u) with distance u from P and the angle at g in T(u) is convex. 


5. Angular metric and structure of non-compact metric surfaces. We now 
assume that an angular measure has been defined for the system of geodesics 
of a G-surface M of finite connectivity. With the notations of the preceding 
section we associate with the points 2; a set of k mutually disjoint tubes 7; 
each bounded by a geodesic polygon P,. 

Let u;>0. By Cohn-Vossen’s Theorem 7; contains a subtube 7;(u,) 
bounded by a geodesic polygon R;(u;) such that R,(u;) is either a closed geo- 
desic, or all angles of R;(u;) measured in 7;(u,;) are concave, or R,(u,;) has 
exactly one convex angle whose measure in 7;(u,) is not zero because R;(u,) is 
a Jordan curve. Let k’ (<k) denote the number of the R,(u,;) with a convex 
angle. 

Call G the compact domain on M bounded by the R,(u,). Since concave 
(convex) angles of R;(u;) measured in T;(u;) are convex (concave) when mea- 
sured in G, the relation (11) yields 
(17) C(G) & 2nX(G)+ k's, 
where the equality sign holds only when all R;(u;) are closed geodesics. 

It is well-known (see [8, pp. 145, 147]) that 
x = {; — (2p +k), p > O if M is orientable, 
2—(p+k), p 2 1 if M is non-orientable, 
where p is the genus of M or G. 

Hence for non-compact M (that is k 2 1) and C(G)2 0 only the following 
cases are possible. If M is orientable, then p = 0 and 1) k = 1, k’'= 0, 1; 
2) k= 2,k’=0,1,2; 3) k= k’=3. If M is not orientable then p = 1, 
k=1,k’=0, 1. 

Taking first only k into account we find in addition to Theorem (13): 

(18) THEOREM. A non-compact G-surface with non-negative curvature is homeo- 
morphic to a plane, a cylinder, a sphere with three holes, or a Moebius strip. 

This agrees again with the known facts regarding Riemann spaces, except 
for the sphere with three holes. It may therefore be of interest to discuss this 
exception in some detail. 

For that and other purposes we divide the tubes, following Cohn-Vossen, 
into three categories. Let T be a tube bounded by the simple closed geodesic 
polygon P, 8 the greatest lower bound of the length of all curves homotopic to 
Pon T. We call minimal sequence a sequence of curves on T homotopic to 
P whose length tends to £. 

If there is no bounded minimal sequence, we call T contracting. 

If no subtube of T is contracting we call T expanding. 

If T is neither contracting nor expanding we call T bulging.* 





*Cohn-Vossen calls a contracting tube a Schaft, and uses Kelch for both bulging and ex- 
panding tubes. The latter are distinguished as “eigentliche Kelche.” 
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The following facts are obvious (Compare [4, §18]): 


(19) A subtube of a contracting tube is contracting. 
(20) A subtube of an expanding tube is expanding. 
(21) A subtube of a bulging tube which is sufficiently far away is contracting. 

An expanding or bulging tube contains a bounded minimal sequence. This 
sequence contains a converging subsequence which tends to a curve R homo- 
topic to P of length 8 (see [1, Sec. I.1]). If the distance of R from P is uw’ then 
R(u)=R for every u > u’. By the Theorem of Cohn-Vossen R(x) is either a 
closed geodesic or all its angles measured in T(u) are concave. 

In the preceding discussion k’ may therefore be interpreted as the number 
of contracting tubes and we see: 

(22) A sphere with three holes and non-negative curvature has only contracting 
tubes and the angle of at least one R;(u;) measured in 7;(u,;) must be less than 
4/3. 

Cohn-Vossen proves that u can be chosen such that the angle of an R(u) 
on a contracting tube is as close to x as desired. This is not true for general 
angular metrics. 

An instructive example can be obtained as follows: In the ordinary space 
consider the surface M of revolution z =(x*+ y*)~'. It is homeomorphic to 
a cylinder or a sphere with two holes, one corresponding to z = @, the other 
toz=0. If P; and P; are two simple closed geodesic polygons associated 
with those holes as in the beginning of this section, say P; toz = © and P, 
to z = 0, and 7; is the tube bounded by P;, then P; is contracting and P, is 
expanding. Well-known facts on geodesics on surfaces of revolution yield 
readily that the R:(u) have all exactly one convex angle D(u) in T,(u) whose 
vertex g(u) has distance u from P. Because M is a surface of revolution and 
the meridians are geodesics the g(u) either lie, or can be assumed to lie, on one 
meridian. 

Let a(u) be the ordinary radian measure of D(u); by Cohn-Vossen’s already 
mentioned result a(u) — x for u—> ©. We now define an angular measure at 
q(u) as follows. If D(u)=(ti, t2)"™ let D be the straight angle of the form 
(ti, t’;) that contains D(u). For any reD let a(r) be the ordinary radian measure 
of (t,, 1)", so that a(r2)= a(u) and define ¢(r) by 


4c) = { da(r) for reD(u), O< 6 <1, 
$a(u)+(a(t)— a(u)) - (# — da(u)) - (x — a(u))™ 
for r<eD — D(u). 
We use ¢(r) as at the end of Sec. 1 to define an angular measure at the point 
q(u). 

For points on the same parallel circle as g(u) we define angular measure in 
an obvious way by rotation of g(u) about the z-axis. On the remainder of M 
we use the ordinary angular metric. Then the new angular metric is con- 
tinuous on M except on the parallel circle corresponding to u—> 0+. It can 
easily be smoothed out there. 

Then D(u) = dx for all u > 0, so that D(u) does not approach x for u— ©. 
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By the same method a sphere with three contracting tubes can be constructed 
for which the angles of all R;(u;) are less than 2/3, so that (22) cannot be 
improved without a new condition on the angular metric. The example shows 
also in which direction such a condition has to go: 

An angular metric is called uniform on a subset G of M if two positive 
functions 5(e) and p(p, «), where 0 < ¢ < 1 and eG, exist, such that the rela- 
tions 0 < aip = dep < p(p, €) and a;a2/(a;p + paz) 2 1 — 4(e) imply for peG 
that |aipa,| > r —«. 

The uniformity is contained in the requirement that 4(e) is independent of p. 
The usual angular metric of a Riemann space is uniform, because |a,pa2| > 
2 arc cos [(1 — 8)/2] for a;— pand a,a2/(a;p + par)= 1 — 8. According to 
Cohn-Vossen (16) may be completed by 


(23) THEoREM. If the angular metric on the tube T is uniform, then for a suitable 
uo> 0 the curve R(uo) is either a closed geodesic, or all angles of R(tuo) measured 
in T(uo) are concave, or the angle of R(uo) at its only vertex q is at least x — «. 

Proof. Consider the function 

f(u) = A(u)+ 28(e)u, u 2 1. 
Since \(u) is non-negative and continuous, f(u) reaches its minimum at some 
value up (2 1). Therefore 
(uot h)+ 25(€)(uo+ hk) > A(uo) + 25(€)uo, for h > 0, 
(23a) 25(e)h > A(u0o)— A(uo+ kh), for h > O. 

If R(uo) is not a closed geodesic or its angles are not concave, let go be the 
vertex of R(uo). If t(go, a*1), t(go, @*2) are proper segments on the legs of the 
angle at go, let (qoa,a;*), i = 1, 2, with h = qoai< p(qo, «). 

Consider the curve R’ originating from R(uo) by replacing t(@:, go) U t(go, a2) 
by t(a1, a2). The distance of R’ from P is at most up+ h. Therefore \(uo+ A) 
{X= length of R’ and 
(24) A(tu0) — A(uto+ h) > Ato) — A’ = Zh — ayae= 2h(1 — asa2/2h), 
and (23a) and (24) yield 


8(e)2 1 — aya2(aigo+ gots); 

hence |agoae| > x — € by the definition of 8(€). 

From (22) and (23) we find 
(25) A sphere with three holes and uniform angular metric cannot have non- 
negative curvature. 

Other well-known theorems can be proved under these general conditions. 
We mention only one example from Hadamard’s theory (see [6)): 
(26) On a G-surface M with -negative curvature a class of freely homotopic 
curves contains at most one closed geodesic. 


The universal covering space of a G-surface M is again a G-surface M’ (see 
[2, Sec. 13]). An angular metric on M induces an angular metric on M’. If 
M has negative curvature, then M’ has negative curvature with respect to 
this induced metric. 
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If M contained two freely homotopic closed geodesics g; and g:, draw a 
segment t from a point #; of g; to a point p, of g2. The figure consisting of 
@1, G2 and t is image of a quadrangle in M’ whose angle sum is 2x. Because 
of the additivity of the excess M’ must contain arbitrarily small non-degenerate 
triangles with non-negative excess, but then M would contain such triangles. 


6. The integral curvature of non-compact surfaces. A polygonal region G 
is the closure of an open set on M whose boundary B (if any) is locally a simple 
geodesic polygon. That means: if p is any point on B then a geodesic triangle 
abc exists which contains in its interior J and such that the intersection of B 
with the closure of J decomposes J and consists of two segments t(p, x), t(p, y). 

For compact G the integral curvature C(G) was defined in Sec. 3. For 
general G we proceed as follows: Let G:C G2C ... be a sequence of compact 
polygonal regions with V G;= G and the further property that a sequence of 
points PieG,;, ,— Gn,, where {n;} is any increasing sequence of positive in- 
tegers, has no accumulation point. If lim C(G,) exists (+ © admitted) it is 
independent of the particular sequence {G,,} and is called the integral curvature 
of G. 

The condition that {P;} has no accumulation point implies for compact G 
that G,= G for large n, so that the present definition of C(G) agrees with the 
previous one. It is necessary to add some such condition because G = UG, 
implies C(G)= lim C(G,) in general only if C(G) can be extended to a com- 
pletely additive set function (compare Sec. 7). 

If M is a G-surface of finite connectivity and the tubes 7; are defined as in 
the beginning of Sec. 5 and H is the compact domain on M bounded by the T;, 
then X(H) is independent of the choice of the 7; and —X(H) is called the 
characteristic —X(M) of M. 

If C(M) exists, it may be evaluated as follows: Let 7," be a sequence of 
subtubes of 7; with T;*C T;"" and (\ T;*= 0. If H® denotes the compact 


n 

domain on M bounded by T;",..., 7,” then 

C(M) = lim C(H"). 
Since any tube, in particular 7;", contains a subtube bounded by a polygon 
R(u) as constructed in (16), it follows from (11) that 
(27) C(M) <& 2eX(M)+ kr, 
provided C(M) exists. The discussion preceding (22) yields 
(28) C(M)2 2xX(M) if M has no expanding tubes. 

An application of (27) is 
(29) A non-compact surface with non-negative, but not identically vanishing curva- 
ture is homeomor phic to a plane. 

For there is a triangle on M with positive excess. This triangle contains 
then a triangle abc with positive excess which is so small that the images of 
abc on the universal covering surface M’ of M are disjoint. If M were nota 
plane M’ would have infinitely many sheets, and in each a copy of abc. The 
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integral curvature of M’, which exists because M’ has non-negative curvature, 
is therefore @. But this contradicts (27). 

Finer results than (27) can be obtained if the angular metric on M is uniform: 

(30) If M has an integral curvature and a uniform angular metric then 
C(M)< 2"X(M). 
The equality sign holds if M possesses no expanding tubes. 

For if the previous notations are used, then 7," carries a polygon R,"(u;) 
with the properties described in (23). If H" is the compact domain bounded 
by the R;"(u,;") then by (11) 

C(H") & 2nX(M)+ k’e. 
The remark about the equality sign follows from (28). 

Theorem (23) yields also 
(31) If the tube T with boundary P has an integral curvature and a uniform ang- 
ular metric then 

C(T)< — 2(e — ai), 
where a; are the angles of P measured in T. 

If C(T) exists, then every subtube of T has an integral curvature. If T is 
contracting or bulging, then it contains T(u) bounded by an R(u) which has 
one vertex q with a convex angle in T or is a closed geodesic. (31) yields then 
C(T(u))< 0. Therefore 
(32) A tube with positive curvature and a uniform angular metric is expanding. 

Cohn-Vossen proves this for tubes with non-negative curvature, but the 
Riemannian character of his metric is essential for this refinement. 

We next prove a theorem which is similar to (31) and is found in the paper 
[5] of Cohn-Vossen. 


@ 
(33) On M let Q be an open Jordan curve of the form  t(ai, Gis1), Gi Gins, 


i=u—@ 

and such that only a finite number of its angles are not straight. Assume, more- 
over, that Q bounds on M a domain G homeomorphic to a halfplane, and that each 
subarc of Q is a shortest connection of its endpoints in G. If a:,..., an are the 
angles at the proper vertices of Q in G and G has a uniform angular metric and an 
integral curvature, then 


c(G)< - , = ed. 


Proof. Let G’ be any simply connected domain in G bounded by a subarc 
Q’ of Q which contains all m vertices of Q and a simple geodesic polygon Q”’ in 
G connecting the two endpoints of Q’. Let peQ’ and let q:(#), g2(t) be the two 
points of Q for which the subarcs from to q;(t) of Q have length ¢. Fora 
proper choice of ¢’ the points g,(t) will lie on Q — Q’ fort 2 t’. Let p(t) bea 
shortest connection of length A(¢) of g:(f) and g2(t) in (G — G’)UQ”. By the 
minimum property of Q 
(34) A(t) 2 2t. 


As in the proof of (16) it is seen that p(¢) is a simple geodesic polygon whose 
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proper vertices, if any, coincide with vertices of Q’’ and such that the corres- 
ponding angles are convex if measured in the domain G(t) bounded by p(t) and 
the subarc from q:(¢) to g2(t) of Q@. By (11) 


C(G(t)) < 2x —(e — Bi(t)—(e — Ba(t))— Z(e — ai), 
where §;(¢) is the angle of p(t) and Q at g; measured in G(#). Due to the arbi- 
trariness of G’ the theorem is proved if a to? ?¢’ exists for which Bi(to) < «. Let 
k = 28(e), where 4(e) is the function entering the definition of a uniform angular 
metric. Then because of (34) 
(35) Mi)— 22+ kt, fort o. 


The triangle inequality implies that \(¢) is continuous, therefore the left side 
of (35) reaches a minimum at some value f92 ¢’. Then 

A(to+ kh) — 2(to+ h)+ k(to+ h)— A(to) + WZo— kto> 0, for h > 0, or 
(36) A(to+ h)— A(t) 2 h(2—k), for h > 0. 


On the other hand if a’;, a’’;, a’s;C p(t), a’ s®C Olte), lie on the legs of the angle 
in (G — G(to))U p(to) at gi(to) (that means \a’ sqi(to)a”’ «| = © — Bi(to)) and 
satisfy the relations 


a’ qi(te) = a” iqi(to) = h < min p(qi(to), €) 


then p(to+ h) is at most as long as the polygon originating from p(to) by re- 
placing t(gi(to), a’;) by t(a’;, a”;). Therefore 
A(to+ h) < A(to)+ (@’1a"1— h) + (a'20""2— h) 
which yields together with (36) 
h(2 — k)< a';a";— h + '2:0":— h, or 
(1 — aa" ,/2h) + (1 — a’2a"’2/2h) { k/2 = &(e). 
Since 1 — a’;a";/2h 2 0 it follows that 1 — a’;a”;/2h < &(e) and from the 
definition of 6(¢) that 
T°=-_ Bi(to) = \a’ qi(to)a’’ s| 2 f<-¢ q.e.d. 
If in addition to the assumptions of (33) Q is a straight line (see (3, p. 232]) 
then there are no corners, hence C(G) < 0. Therefore 
(37) A plane with positive curvature does not contain a straight line. 
If the assumption that every subarc of Q is a shortest connection in G is 
omitted, Cohn-Vossen proves that 
(38) C(G) & « — U(x — aj). 


In general spaces the inequality C(G) < 2x — 2(x — a,) is trivial, but the 
refinement from 2 to x rests on the fact that in Riemannian geometry perpen- 
dicular directions form the angle +/2. This fact has no analogue in general 
Finsler spaces, no matter how the angular metric is defined, because perpen- 
dicularity is not symmetric. Consequently there is no reason to believe that 
(38) holds with a suitable definition of angular measure, unless perpendicularity 
is symmetric, although the author did not try to construct an example because 
this would obviously be very laborious. 
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We conclude the analysis of the validity of Riemannian methods in general 
spaces, which could be continued almost ad libitum, by mentioning that the 
proofs of the following two interesting results of Cohn-Vossen [5] hold without 
any change: 

Let M be a plane with positive curvature and uniform angular metric. Then 
every point of M lies on at least one geodesic without multiple points. If a geodesic 
g has multiple points, then it contains exactly one 1-gon P, moreover g — P lies 
in the exterior of P and consists of two branches without multiple points (but the 
two branches may intersect each other). 


7. The integral curvature as set function. On a surface M with a system S 
of geodesics and an angular metric as defined in Sec. 3, let Fo be the collection 
of the following sets: the empty set, the points, the segments without end- 
points (l-cells), the interiors of the non-degenerate triangles (2-cells). The 
excess is called completely additive if for any representation of a 2-cell abc as 
union of a countable number of disjoint 2-cells a,b,c, and points and 1-cells on 
the boundaries of the a,b,c, 


e(abc) = > e(arb»cr). 


The unions ¢o of a finite number of disjoint elements in Fy form a field F;. 
If a,b,c, are the two cells of a given set ceF; we put 


C(o) = Ze(a,b,c,). 


If C(c) is bounded on every bounded subset of M and the excess is completely 
additive, then C(c) can be extended to a completely additive set function on the 
o-field F of all Borel setson M. Moreover the extended set function is bounded 
on every bounded subset of M with the same bounds as the old function.® 

Let a measure m be defined on M for which segments have measure 0, and 
such that every bounded measurable set on M has finite measure. We are 
going to prove the theorem 
(39) If for every bounded set B on M a number 8(B) exists such that for any 2-cell 
abc in B 

|«(abc)| < 8(B)m(abc) 
then ¢ is completely additive, C(a) is bounded on every bounded subset of M and 
absolutely continuous on F with respect to m. 


Proof. Leto bea set in F; which lies in the given set B and a,b,c, the 2-cells 
of «. Then 


C(o) < Z\e(a,b,c,)| < (B) - Zm(a,b,c,) < B(B)m(c) < B(B)m(B), 


so that C(c) is bounded in B for ceF;. 
If abc is the union of the disjoint 2-cells a,b,c,, vy = 1, 2,... and points and 


*The arguments which lead to these conclusions are implicitly contained in many modern 
treatments of set functions. For those who are able to read Danish an unusually clear expo- 
sition is available in Jessen (7, part 3] which also determined the present formulation. 
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l-cells 5; on the boundaries of the a,b,c,, then m(abc) = Em(a,b,c,) because 
m(é;) = 0. For a given e > 0 we can therefore find an n(e) such that 
n(e) 
m(abc) — > m(a,b,c,) < €/B(abc). 
v=l 
n(e) 
Then abc — > a,b,c, is the sum of a finite number of 2-cells a’;b’;c’;, 


y=l 
i= 1,..., mand a finite number of points and l1-cells. Therefore 
ad m n(e) 
X e(a’sb’ ic’) < B(abc) | m(a’ xb’ ;c’;) = B(abc)m(abe — ¥ a,b,c,) X « 
i=l i=1 r=l 


which shows that e(abc) = > ¢(a,b,c,) or that ¢ is completely additive. 


By the preceding remarks and the first part of this proof C(c), ceF, can be 
extended to a completely additive function on F with the same bounds. It 
then follows that 

|C(c)| < 8(B)m(c) for oeF and o C B. 
Therefore m(c) = 0 implies C(¢) = 0 so that C(c) is absolutely continuous. 

Under the hypotheses of the theorem, C(c) is therefore the indefinite integral 
f. f(p) of a function f(p) with respect to the measure m. This does not yet 
assign a definite value to f(p) at any given point since f(p) can be changed at 
will in a set of measure 0. This indefiniteness can be eliminated if sufficient 
restrictions on the angular metric and on the measure guarantee that there 
is at least (and then exactly) one continuous f(p). But it seems more worth- 
while to discuss these questions in connection with a specific angular measure 
in a Finsler space. 
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THE DENSITY OF REDUCIBLE INTEGERS 
S. D. CHOWLA AND JOHN TODD 


Introduction. The concept of a reducible integer was introduced recently 
[3] : if P(m) denotes the greatest prime factor of m then n is said to be reducible 
if P(1 + m*)< 2n. The reason for the term is that reducibility is a condition 
necessary and sufficient for the existence of a relation of the form 


r 


arctan n = 2 fi arctan n; 
is 
where the f; are integers and the m; positive integers less than n. J. C. P. 
Miller pointed out to us the regularity of the distribution of the reducible 
integers (less than 600). In collaboration with Dr. J. W. Wrench, using his 
tables of factors of 1 + n’, we carried the count still further, and observed the 
same regularity. The following conjecture suggested itself: 
ol “‘Reducible integers have a density about 0.3.” 
We have not been able to make very much headway with this but have 
succeeded in establishing the following: 

TuHEoreM A. The density of the set of integers n for which P(n)< 2n' is 
1 — log 2 =.3069.... 

This note contains a proof of this theorem, and a table summarizing the 
numerical evidence in support of C. 

1. Numerical evidence. We give here a summary of the numerical evi- 
dence relating to the conjecture C together with corresponding results related 
to Theorem A. The table below gives, in each range (1 + 100 m, 100(m + 1)), 
for nm = 0(1)49, on the right, the number of reducible integers in that range, 
and on the left, the number of integers in that range which satisfy P(n) < 2n'. 

Totals in the various chiliads and a grand total for the complete range 
(1-5000) are given in the last line of the table. 


0 1000 2000 3000 4000 

1-100 (29,57) (31,43) (29,43) (83,41) (29, 42) 
101-200 (29,50) (25,43) (30,42) (28,43) (28, 40) 
201-300 (28,47) (33,44) (23,42) (23,43) (27,41) 
301-400 (26,45) (28,41) (32,41) (82,43) (31, 40) 
401-500 (30,45) (31,44) (28,44) (29,38) (27, 42) 
501-600 (30,44) (23,44) (32,39) (32,41) (38, 39) 
601-700 (30,44) (27,40) (26,43) (25,40) (30,41) 
701-800 (29,44) (34,43) (32,41) (30,43) (35,39) 
801-900 (27,44) (28,45) (27,42) (29,40) (30, 43) 
901-1000 (23,42) (31,39) (29,41) (19,41) (38, 41) 

(281, 462) (291,426) (288,418) (280,413) (313, 408) (1453, 2127) 


Received April 5, 1948. This work has been supported by the Office of Naval Research 
of the United States Navy Department. 
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2. Proof of Theorem A. It is more convenient to show that the density 
of the integers m for which P(n)> 2n',islog 2. That is, we shall show that 


Ox)= LY 1l~xlog?2; 


P(n)22n4 
to do this we establish the two following results: 
Ay. Qix)= LD 1~x log 2; 
reanat 
A Q2(x) = Q(x) — Q(x) = >» 1 = o(x). 


and < P(n) <2st 

2.1. Proof of Ai. This is carried out by a modification of a method used 
recently [1] to evaluate lim x'R,(x) where R,(x) is the number of integers 
n < x for which P(n)> x*. 

For any p the number of integers m < x which are multiples of p is [x/p). 
In Q:(x) we consider only primes p = P(n)> 2x*: for such primes the residual 
factor (n/p) < 4 x*< p and so every multiple of p which does not exceed x 
has p for its greatest prime factor. Hence 


Qi(x)= 2 [x/p 


223<pSxz 
= > {(«/p) + O(1)} 
axb< psx 
=x 2 p+ O(x/log x), 
axb<p<x 
since SY 1< EY 1 = O(x/log x). 
axb<p<x otz 
It is, however, well known [2, pp. 100-102] that 
B. > p= log log x — 1 + O(1/log x) 
osx 


where / is a certain constant. Hence 

x10,(x) log log x — log log 2x'+ o0(1) 
log {(log x)/(} log x + log 2)} + o(1) 
log 2 + o(1), 


which establishes A. 


2.2. Proof of Az. This is carried out in the following manner. First, it 
will be sufficient to restrict the values of m considered to the range 
x/(log x)*Sn< x, 
for this implies a change in the sum of O(x/(log x)*) = o(x). Secondly, we do 
not decrease the sum if we replace 2n', the variable limit in the lower inequality, 
by its smallest value 2x*/log x. Thirdly, we do not decrease the sum by now 
allowing m to cover the full range 1 < » <x. Thus it will be sufficient to 


show that Qs(x) -_ 2 l= o(x). 


au >= 
(224 /1og x) < p(n) < 2x4 
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In order that an integer should contribute to Q, it is necessary that it should 
have a prime factor p in the range (2x*/log x, 2x"). For p fixed the number 
of such n is [x/p]. Hence 
Q< y [x/p] . 
(2x4/log x)<p<2xt 
(It is possible for an integer » < x to have two factors in the range and so we 
must allow for inequality, which was not so in the case of Q,.) 
We now proceed as before: 
Q:(x) < ) [x/p] = x z p+ O(x*/log x) 
(2x4/log x) < p< 2x4 (2x4/log x) < p< 2x4 
= x {log log 2x*— log log (2x*/log x)} + O(x/log x), 
using B. Since 
log log 2x'— log log (2x'/log x) 
= log {(4 log x + log 2)/(} log x + log 2 — log log x)} 
log [{1 +(log 4)/(log x)} {1 +(log 4 — 2 log log x)/log x} ] 
log {1 + O(1/log x))(1 + O(log log x/log x)} 
O(log log x/log x) = o(1), 
the proof of A; is complete. 


3. Possible generalizations. It is clear that 2n* in Theorem A can be 
replaced by An’ for any A > 1 without affecting the conclusion. 

Similar arguments show that the density of the integers m for which P(n) > 
An* (§ <a <1,A > 1) is exactly log a. 

The case when a < } requires more careful study along the lines indicated 
in [1] and it can be shown that the device used here (replacing a summation 
over 1 < m < x by one over (x/(log x)*< nm < x) will enable the density to 
be evaluated explicitly in this case, too. 

It is clear that an estimate for the error term 

x"O(x) — log 2 
is 
O(log log x/log x), 
and this explains the slowness of the convergence apparent in the table. 
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ON A THEOREM OF LATIMER AND MACDUFFEE 


OLGA TAUSSKY 


THE matrix solutions of an irreducible algebraic equation with integral co- 
efficients were studied by Latimer and MacDuffee.' They considered matrices 
with rational integers as elements. If A is such a matrix, then all matrices 
of the “class” S~ A S will again be solutions if S is a matrix of determinant + 1. 
On the other hand, in generai all solutions cannot be derived in this way from 
one solution only. It was in fact shown that the number of classes of matrix 
solutions coincides with the number of different classes of ideals in the ring 
generated by an algebraic root of the same equation. Although this result 
is of interest in many different branches of mathematics it is not generally 
known. It seems particularly often required for periodic matrices.” 

Latimer and MacDuffee actually dealt with the more general case when the 
equation was reducible. By restriction to irreducible equations only, a very 
simple proof can be obtained. 

In what follows f(x) =0 is an irreducible algebraic equation of degree m with 
integral coefficients, a one of its algebraic roots, A =(a4) an m X nm matrix 
with rational integers as elements which satisfies f(x)= 0 and S is a matrix 
with rational integers as elements and determinant + 1. 


THEOREM 1. The algebraic number a is a characteristic root of the matrix A 
and the components of the corresponding characteristic vector (a,..., an) can be 
chosen to form the basis of an ideal in the ring formed by the polynomials in a 
with rational integers as coefficients. 


Proof. Since f(x) is assumed irreducible it follows that it is the character- 
istic and the minimum polynomial of A and that a is a characteristic root of A. 
Since in this case the characteristic roots of A are ali simple, the corresponding 
characteristic vector is uniquely determined apart from a factor of proportion- 
ality. Since 
(1) a(ai,.--,@n) = A(as,..., Gn) 


we may take for a; the cofactor of the ith element in a fixed row of the deter- 


Received July 20, 1948. 
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minant |ai,— ad%|. This is a polynomial in a with rational integral coeffi- 
cients. From (1) it follows that 
a‘(a1,...,@n) = A*(a;,..., Gn), (@=0,...,”—1). 
Since the numbers 1, a,..., a"~' form a basis for the ring in question, it is 
proved that the set of numbers 
Qiait. . .+ Anan 
where a; are rational integers, forms an ideal. 


THEOREM 2. Two ideals determined (as in Theorem 1) from the same matrix 
A belong to the same ideal class. 

Proof. Since the elements of the basis of an ideal in Theorem 1 are uniquely 
determined apart from a common multiplier it follows that any two such ideals 
must be equivalent—as usual two ideals a and 6 are said to be equivalent 
or belong to the same class if two elements, a, 8 in the ring exist such that 

aa = bf. 

THEOREM 3. To every ideal (w:,...,@n) im the ring generated by a there 


corresponds a matrix X with rational integers as elements which satisfies f(X) = 0 
and such that 


a(wi,..-,@n) = X(w1,..., Wa). 
Proof. Since (w1,..., @n,) is an ideal there must exist a relation 
(2) a(wi,..-,@n) = X(w,..., Wa) 
where X is a matrix with rational integral elements. From (2) follows 
a*(wi,...,@n) = X*(wi,..., Wn), (@=0,...,”—1). 


This implies 
f(a)(ws,..., @n) = f(X)(@1,..., @n) = 0. 
Since f(X) is also a matrix with rational elements and since the relations 
f(X) (or, . . ., an) = 0 
hold in the fields generated by the conjugate roots of a and |w;”| 0 it follows 
that 
f(X) = 0. 

THEOREM 4. The matrix X in Theorem 3 is uniquely determined apart from 

a transformation SX S“. 


Proof. Ifa different basis for the ideal were chosen it would be of the form 
S(wi,..., @n) with |S| = +1. There would then be a relation 


aS(wi,..-,@n) = YVS(w1,..., Wa). 
On the other hand, in virtue of (2) 
aS(wi,.-.-.,@n) = SX(w1,... 5 Wa)- 
Hence by the argument used at the end of the proof of Theorem 3: 
SX = YS or 
Y = SXS". 
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The Theorems 1—4 show that there is a 1-1 correspondence between the classes 
of matrices and the ideal classes. 

It may be pointed out that the matrices S for which SAS“'= A play a role 
similar to the units in the algebraic number fields. Such a matrix is in fact 
a polynomial in A, and since its determinant is + 1 it is a unit in the field 
generated by A. 
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See e.g. J. H. M. Wedderburn, “Lectures on Matrices,” Amer. Math. Soc. Colloquium 
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CONGRUENCE RELATIONS BETWEEN THE TRACES 
OF MATRIX POWERS 


J. S. FRAME 


1. Introduction. Let A be a matrix of finite order m and finite degree d, 
whose characteristic roots are certain n roots of unity a;, a2... , ag. We wish 
to prove a congruence (6) between the traces (tr) of certain powers of A, which 
is suggested by two somewhat simpler congruences (1) and (3). 


First, if tr (A) is a rational integer, it is easy to establish the familiar con- 
gruence 


(1) tr(A) = tr(A”) (mod p), p prime, 
even though tr(A”) may not itself be rational. For we have 
d » 6 
(2) [tr(A)]? = | = x | = 2 a,?+ p(...) = tr(A”) + p(...) 
where (...) denotes an algebraic integer. The left-hand members of (1) and 


(2) are rational integers which are congruent (mod p) by Fermat's theorem. 


The right-hand members are explicitly congruent (mod p). Hence (1) follows 
from (2). 


Secondly, for any integer a, we have 


(3) a” = a? (mod p*), if # isa prime power > 1. 
Equation (3) is trivial if a is divisible by p. Otherwise it can be established 
easily by setting m = p* in the well-known Euler congruence 
(4) a*™ = 1 (mod m), for (a, m) = 1, 
where ¢(m) is the Euler ¢-function, and ¢(p*) = p*— p**". 
It is our purpose to prove a congruence relation (mod p*), which generalizes 


(1) and is similar to (3), between the traces of certain powers of a matrix A 


of finite order—or, in other words, between certain sums of powers of roots 
of unity. 


THEOREM. Let S(m) denote the trace of the m power of a matrix A of finite 
order n and finite degree S(O), and assume that A is such that 


(5) S(k) = S(1), for all k such that (k,n) = 1. 
Then 
(6) S(p?) = S(p’*) (mod 7’). 


We note that condition (5) implies that A has a rational integral trace, but 
that not every matrix with rational integral trace satisfies (5). 


Received July 21, 1948. 
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2. Proof of the theorem.' Let us define a ‘“‘p*-set” to be a set of roots of 
unity such that the sum of its #**" powers are congruent to the sum of its p*~"** 
powers mod #* as in (6). We note that the negative of any root of unity is 
also a root of unity. 


Lemma 1. The set of all the n distinct n“ roots of unity is a p*-set. 


Proof. Denoting the sum of m‘" powers by S,(m) we have 
(7) S.(p*) = n, if n divides p’, 
0, if m does not divide p*. 
Hence 
(8) S.(p*) S,(p°~") pf ifn = t, 


= 0 otherwise. 


Lemma 2. If one p*-set is included as a subset of a larger p*-set, the 
difference of the two p*-sets is also a p*-set. Furthermore, any set of roots of 
unity which is made up of two or more p*-sets is also a p*-set. 


Proof. If each of two or more quantities S(p*) — S(p*~') is congruent to 0, 
so is their sum or difference. 


Lemma 3. The set of ¢(m) primitive n‘® roots of unity is a p*-set. 


For prime n the lemma is a special case of Lemma 2. Assuming as induction 
hypothesis that the lemma is true for all » with a smaller number of prime fac- 
tors than n, we show that it is also true for » by applying Lemmas 1 and 2, and 
eliminating from the complete set of  n“” roots all sets of primitive »“” roots 
for each v which is a proper divisor of nm. Only the primitive n‘* roots remain. 
They form a p*-set. 

We observe that condition (5) implies that for any factor yu of m the primitive 
py roots occur as roots of the matrix A with equal multiplicity. Hence by 


Lemmas 2 and 3 the roots of A are a p*-set, so the theorem is established. 


3. Applications of the theorem. In constructing the table of characters 
for a finite group, our theorem may be applied to determine many of the entries. 
For example, the symmetric group of degree 5 and order 120 has irreducible 
representations of degrees 1, 1, 4, 4, 5, 5,6. There are 15 conjugate elements 
of order 2 which are squares of elements of order 4. Hence their traces form 
a vector of unitary squared length 120/15 which is unitary orthogonal to the 
vector (1, 1, 4, 4, 5, 5, 6) and congruent to it (mod 4). The only integral solu- 
tion is (1, 1,0,0,1,1, —2). Similarly for the traces of the 24 elements of order 5 
we have the vector (1, 1, —1, —1, 0,0, 1) as is known by the ordinary modular 
theory (mod 5). Given the numbers of elements in the classes of conjugates, 
the table is completely determined by these congruences. 
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1] am indebted to Professor R. Brauer for some suggestions for shortening my original proof. 





